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is  important  to  the  design  process  to 
forecast  the  degree  of  wave  damping  for  a 
berm  of  given  dimensions  and  mud 
rheology. 

RESEARCH:  The  mud/wave  domain  and 
fluid  mud  were  defined  and  the  equations 
identified  that  describe  generation  and 
erosion  of  fluid  mud.  Water  waves  over 
mud  and  water-column  dynamics  were 
modeled.  Model  simulations  were 
compared  with  laboratory  and  field  data. 


SUMMARY:  As  its  title  indicates,  the 
report  presents  a  short  review  of  iterative 
procedures  and  simple  modeling  approaches 
for  studying  the  basics  of  the  interaction 
between  progressive  water  waves  and  a 
compliant  mud  bed  in  the  shallow  coa.sta! 
environment.  Laboratory  and  field  data  and 
model  simulations  are  discussed;  an 
appendix  briefly  summarizes  additional 
modeling  studies. 
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Summary 


The  basics  of  the  interaction  between  progressive  water  waves  and  a 
compliant  mud  bed  in  the  shallow  coastal  environment  have  been  briefly 
reviewed.  A  key  feature  of  the  interactive  process  is  the  ability  of  waves  to 
fluidi2e  bottom  mud  and  sustain  it  in  that  state  as  long  as  wave  action  con¬ 
tinues.  Compliant  or  fluid  mud  is  a  highly  viscous  medium  that  can  oscillate 
as  waves  pass  over  and  cause  v.ave  heights  to  attenuate  significantly.  At  the 
same  time,  given  a  high  enough  fluid  stress  at  the  mud-water  interface,  panic¬ 
ulate  entrainment  can  occur,  thus  increasing  water  column  turbidity. 

The  wave-mud  interaction  problem,  involving  the  prediction  of  surface 
wave  attenuation,  bottom  mud  motion,  and  interfacial  entrainment  or  erosion, 
is  simply  treated  here  as  one  concerned  primarily  with  vertical  exchanges  of 
momentum  and  sediment  mass.  A  dichotomy  inherent  in  such  a  treatment 
arises  from  the  need  to  assume  mud  to  be  a  continuum  when  simulating  wave 
attenuation  and  mud  motion,  while  considering  the  vertical  transport  of  sedi¬ 
ment  in  the  water  column  to  be  a  two-phased  problem  amenable  to  classical 
approaches  in  sediment  transport. 

At  the  center  of  the  continuum  approach  is  the  requirement  to  describe  the 
constitutive  relations  characterizing  mud  rheology.  Depending  on  the  prop¬ 
erties  of  the  constituent  sediment  and  the  ambient  fluid,  mud  in  general  can 
range  from  being  a  highly  rigid  and  weakly  viscous  material  to  one  that  can 
be  approximated  as  a  purely  viscous  fluid.  In  that  context,  fluid  mud  is  better 
defined  as  fluidlike  mr'^,  since  it  is  not  always  wholly  devoid  of  rigidity. 

As  wave  action  proceeds,  the  viscoelastic  (or  poroelastic)  properties  of 
mud  tend  to  change  with  time  over  scales  that  can  be  two  to  three  orders 
longer  than  the  typical  wave  period.  The  resulting  feedbacks  in  terms  of  the 
influence  of  this  change  on  wave  propagation  are  not  always  easy  to  quantify. 
Yet,  over  the  past  three  to  four  decades  simple  hydrodynamic  models  have 
been  developed  to  simulate  wave  attenuation,  mud  motion,  and  interfacial 
erosion  in  the  prototype  environment  with  an  acceptable  degree  of  accuracy. 

Conventionally,  model  applications  have  dealt  mainly  with  investigations  of 
natural  phenomena  including  predictions  of  coastal  wave  heights  and  turbidity- 
generation  for  understanding  beach  erosion  and  flooding,  planning  of  port  and 
harbor  facilities  including  safe  navigation  channels,  <ind  design  of  offshore 
structures  where  bottom  stability  may  be  a  problem.  In  recent  years 


applications  have  been  extended  to  problems  of  water  quality  and  benthic 
biota,  since  both  can  be  influenced  by  mud  motion  and  associated  constituent 
fluxes.  Most  recently,  engineering  interest  has  grown  in  beneficially  using 
dredged  fine-grained  material  from  navigation  channels  to  create  offshore 
underwater  berms  that  can  absorb  wave  energy  and  thus  act  as  buffers  against 
wave  attack  in  areas  in  the  lee  of  the  berms.  To  that  end,  applications  of 
wave-mud  interaction  modeling  have  been  carried  out  to  simulate  the  degree 
of  wave  damping  that  will  occur  for  a  berm  of  given  dimensions  and  mud 
rheology.  Thus,  in  turn  this  type  of  modeling  can  be  potentially  useful  in 
developing  guidelines  for  designing  wave  energy  absorbing  berms  in  future. 


INTRODUCTION 


1 . 1  MUD-WAVE  DYNAMICS  IN  ENGINEERING  AND  SCIENTIFIC  APPLICATIONS 
Advancements  in  numerous  scientific  and  engineering  applications  in  the  shallow  marine  environment  are 

contingent  upon  a  better  understanding  of  the  physical  processes  and  associated  feedback  mechanisms  that  are 
inherent  in  the  dynamical  response  of  a  muddy  bottom  subjected  to  wave  action.  Some  relevant  areas  and  their 
descriptions  are  briefly  noted  in  Table  1.  They  highlight  the  wide  ranging  scientific  and  engineering  concerns  in 
which  mud  dynamics  plays  a  key  role,  including  coastal  fisheries  and  agriculture,  shore  and  channel  protection,  oil 
and  gas  exploration,  water  quality,  shallow  marine  habitats  and  so  on.  Despite  the  diversity  and  complexity  of  the 
physical,  physicochemical  and  biological  processes  involved  in  these  situations,  in  most  cases  the  processes  and 
feedbacks  related  to  the  physics  of  mud-wave  interaction  can  be  identified  in  a  simplified  manner,  and  is  the  focus 
of  this  report. 

1.2  PROCESSES  AND  FEEDBACKS 

The  response  of  a  mud  bed  to  wave  forcing  is  schematized  in  Fig.  1.  A  brief  description  of  the  numbered 
pathways  is  as  follows;  1,2)  water  waves  determine  the  flow  field;  3,4)  flow  field  and  bed  properties  together 
govern  the  character  and  the  dynamics  of  the  interface;  5,6)  wave  loading,  consolidation,  fluidization  and  thixotropy 


Figure  I.  Schematic  description  of  the  processes  and  feedbacks  related  to  the  response  of  mud  bed  to 
water  waves.  Box  with  dashed  line  represents  the  mud-water  system,  inverted  triangle  represents 
forcing,  rectangles  are  components  of  the  mud-water  system,  ellipses  are  process  transfer  functions  and 
circles  represent  manifested  responses  (adapted  from  Maa  and  Mehta,  1989). 
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erosion;  11,12)  shear  (and  normal)  stresses  together  with  mud  properties  determine  mud  motion;  13,14)  mud  change 
mud  properties  with  time;  7,8)  flow  field  and  interfacial  character  determine  the  inierfacial  shear  stress;  9,10) 
interfacial  shear  stress  and  interfacial  properties  determine  the  rate  of  particle  entrainment,  or  interfacial  properties 
largely  determine  the  rate  of  surface  wave  damping  or  attenuation. 

In  this  report  we  will  examine  the  inter-relationship  between  forcing  and  responses  identified  in  Fig,  1 ,  and 
in  so  doing  inherently  highlight  significant  feedbacks  represented  by  some  of  the  numbered  pathways.  Feedbacks 
shown  tend  to  be  time-dependent,  typically  varying  the  behavior  of  the  system  gradually  in  comparison  with  the 
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forcing  period,  as  the  mud  properties  change  under  continued  wave  action.  Thus,  predictive  approaches  for 
determining  wave  attenuation,  mud  motion  and  erosion  are  dependent  on  a  knowledge  of  the  constitutive  properties 
of  mud,  which  themselves  vary  widely  depending  on  the  mineralogical  composition  of  mud,  and  modulation  of 
inorganic  properties  by  biochemical  and  biophysical  influences.  Consequently,  a  few  important  tests  for 
characterizing  mud  properties  and  transport,  e.g.  erodibility,  must  accompany  any  predictive  effort.  However, 
notwithstanding  site-specificity  that  is  inherent  in  every  problem,  an  attempt  has  been  made  in  the  report  to 
emphasize  those  physical  principles  that  are  common  to  most  problems. 

The  degree  of  importance  of  a  particular  response  of  the  mud-water  system  evidently  depends  on  the  problem 
area  and  the  application  sought.  Thus,  for  example,  wave  attenuation  is  of  primary  significance  for  mitigation  of 
shoreline  erosion  and  coastal  flooding.  Mud  motion  determines  bottom  stability,  hence  the  integrity  of  structural 
foundation.  It  also  controls  the  intake  and  release  of  nutrients  and  contaminants  across  the  mud-water  interface. 
Finally,  coastal  and  estuarine  turbidity  due  to  wave  action  is  the  result  of  particulate  entrainment  at  the  mud /water 
interface.  In  any  event,  as  a  first  step  it  is  essential  to  introduce  basic  definitions  related  to  the  physical  state  of  mud 
in  the  wave  field.  This  aspect  is  considered  next. 
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MUD-WAVE  DOMAIN 


2.1  DEFINITIONS 

Three  previous  definitions  of  mud  are  given  in  Table  2,  two  ocean  science  related  and  the  third  from  coastal 
engineering.  From  the  present  standpoint,  the  coastal  engineering  description  is  perhaps  more  relevant  than  the  other 
two  because  it  refers  to  the  state  of  mud  (fluid-to-plastic)  as  opposed  to  composition  alone.  In  that  context  it  is 
noteworthy  that  wave  attenuation  and  mud  motion  are  both  dependent  on  the  constitutive  relations  characterizing 
mud  rheology,  which  in  turn  depends  on  the  physical  state  of  mud.  On  the  other  hand,  particle  entrainment  at  the 
interface  largely  depends  on  the  structure  of  the  particulate  aggregate  network,  which  is  specified  by  sediment 
composition  for  a  given  composition  of  the  fluid  {ionic  species  in  water,  their  concentration,  and  pH).  Mud  that  is 
predominantly  composed  of  particles  smaller  than  about  20  ftm  is  much  more  cohesive  than  one  in  which  the 
particles  are  coarser  (Mehta  and  Lee,  1994).  Where  it  is  essential  to  make  this  distinction,  we  will  refer  to  the 
former  as  cohesive  mud. 


Table  2.  Selected  definitions  of  mud 


Source 

Field 

Definition 

Hunt  and  Groves  (1965) 

Ocean  science 

Pelagic  or  terrigenous  detrital  material  consisting  of 
particles  smaller  than  sand,  i.e.,  an  undifferentiated 
sediment  made  up  of  particles  mostly  within  the 
silt-clay  range  smaller  than  0.0625  mm. 

Allen  (1972) 

Coastal  engineering 

A  fluid-to-plastic  mixture  of  finely  divided  particles 
of  solid  material  and  water. 

Tver  (1979) 

Ocean  science 

Pelagic  or  terrigenous  detrital  material  consisting 
mostly  of  silt  and  clay-sized  particles  (less  than 

0.06  mm),  but  often  containing  varying  amounts  of 
sand  and/or  organic  materials.  It  is  a  general  term 
applied  to  any  sticky  fine-grained  sediment  whose 
exact  size  classification  has  not  been  determined. 

The  forcing  frequency  (f)  of  conunon  interest  relative  to  gravity  waves  ranges  widely,  from  about  10  Hz  to 
10  ’  Hz.  The  former  corresponds  to  the  transition  whereby  surface  tension  becomes  increasingly  important  as  a 
restoring  force  with  increasing  frequency,  while  the  latter  is  representative  of  the  frequency  of  the  astronomical  tide. 
It  should  be  added,  however,  that  frequencies  that  are  lower  than  tidal,  especially  those  corresponding  to  sub-tidal 
oscillations  having  periods  on  the  order  of  days,  can  also  be  important  in  some  situations,  for  example  with  regard 
to  long  term  erosion  or  accretion  of  fine-grained  sediments  in  wind-forced  bays. 

Given  water  depth,  h,  and  particle  settling  velocity,  w,,  two  characteristic,  frequency  dependent  numbers  that 
together  characterize  the  primary  nature  of  wave  forcing  and  bottom  response  are  hw’/g  and  hw/w^,  where  u=27rf 
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and  g  is  the  acceleration  due  to  gravity.  When  the  first  number  is  less  than  0.3,  fluid  pressure  is  practically 
hydrostatic,  and  the  celerity  of  the  shallow  water  wave,  f=(gh)''^.  As  hw-/g  exceeds  0.3,  dynamic  pressure  effects 
become  increasingly  important  with  increasing  water  depth,  or  frequency.  The  second  number  scales  the  water  depth 
in  relation  to  the  settling  velocity  for  a  given  frequency.  For  fine-grained  sediment,  selecting  w,  =  lO"*  m/s,  h  =  10 
m  and  w=l  rad/s  as  characteristic  values  yields  ha)/Ws=I0^.  which  is  up  to  two  orders  of  magnitude  greater  than 
the  corresponding  value  for  a  sandy  bed,  for  example.  A  manifestation  of  this  difference  is  related  to  the  vertical 
structure  of  density,  which  assumes  considerable  importance  when  the  sediment  is  fine-grained.  Sediment-induced 
stratification  characteristically  occurs  in  this  case  irrespective  of  the  value  of  hwVg,  as  long  as  the  wave  is  not  in 
deep  water  (hw/f>Tr),  since  in  the  latter  case  wave-bottom  mud  interaction  ceases. 

Other  parametric  approaches  have  also  been  proposed  for  characterizing  the  wave-mud  environment  in 
contrast  to  wave-sand  environment.  For  example,  McCave  (1971)  has  defined  a  "wave  effectiveness  parameter"  as 
the  product  of  the  theoretical  instantaneous  sediment  transport  rate  times  wave  frequency.  In  the  southern  North  Sea 
area  examined,  relatively  low  values  of  this  parameter  correlated  with  the  occurrence  of  mud. 

2,2  FLUID  MUD 

A  simple  description  of  the  mud-wave  system  is  given  in  Fig.  2,  in  which  u„,  is  the  amplitude  of  the  horizontal 
wave  orbital  velocity.  The  (bulk)  density  of  the  water-mud  mixture,  p,  varies  from  p,  at  the  surface  to  p^  at  the 
bottom.  The  quantity  is  water  density.  A  sharp  density  gradient,  or  lutocline  (Parker  and  Kirby,  1982),  separates 
the  upper  column  suspension  from  fluidized,  compliant  mud  below.  At  the  bottom  of  the  mobile  suspension  the 
density  rises  from  p^  to  p„  which  is  the  range  over  which  fluid  mud  is  considered  to  occur.  The  bed  below  fluid 
mud  can  undergo  deformation.  This  deformation  in  fact  may  eventually  break  the  inter-particle  or  inter-aggregate 
bonds,  and  thus  change  the  bed,  possessing  a  structured  matrix,  to  fluid  mud.  Below  the  level  at  which  the  depth 
of  penetration  of  the  wave  orbit  practically  ends  the  bed  remains  uninfluenced  by  wave  motion. 

Fluid  mud  characterizing  densities,  p^,  and  p„  are  operationally  defined  as  follows.  The  flux  of  sediment 
settling  within  the  mobile  suspension,  F^^w^p^tfi,  where  4)=ip-pJ/{prPw)  is  the  solids  volume  fraction  (  =  1-n,  where 
n  is  the  porosity)  and  p,  =  granular  density  of  the  sediment.  In  Fig.  3a,  data  from  a  laboratory  settling  column  using 
sediment  from  the  Severn  Estuary  in  the  United  Kingdom  have  been  used  to  plot  F,  against  the  excess  density, 
Ap=p-p„  (adapted  from  Odd  and  Rodger,  1986).  The  peak  value  of  F5  =  F5„,=40  g/m-  s  corresponds  to  the  onset 
of  hindered  settling,  under  which  F,  decreases  with  increasing  Ap  (Mehta,  1990).  Hindered  settling  begins  at 
Ap„=p„,-p*=9  kg/m’,  where  p^,  is  the  value  of  p  when  F5=F5„,.  We  further  note  that  at  the  peak  flux,  dFJdAp-Q. 
Furthermore,  given  z  as  the  vertical  coordinate,  dFJdAp=(d¥Jdz)l(dApldz).  Thus,  since  at  the  lutocline  dApIdz 
tends  to  infinity  the  condition,  3F,/dAp=0,  is  satisfied  there  (horizontal  line  in  Fig.  2),  irrespective  of  the  value  of 
dFJdz.  Thus,  in  reality  Ap„  practically  coincides  with  Ap„=p„-p,..  In  Fig.  3b  the  correspondence  between  the 
lutocline  elevation,  z„,  and  Ap„  is  shown  for  a  measured  density  profile  from  the  Severn.  Note  that  the  .sediment 
remains  suspended  in  this  estuary  largely  as  a  result  of  tidal  flow,  as  opposed  to  wave  action.  Where  waves 
predominate,  the  essential  description  should  remain  unchmged,  however. 
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Figure  2.  Schematic  diagram  showing  the  vertical  variation  of  the  density  of  mud-water  mixture  in  the  wave 
environment.  MWL  =  mean  water  level. 
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Fig.  3.  a)  Settling  flux  variation  with  excess  density  for  sediment  from  the  Severn  Estuary,  United  Kingdom. 
Laboratory  column  data  reported  by  Odd  and  Rodger  ( 1 986);  b)  An  instantaneous  excess  density  profile  showing 
a  marked  lutocline  in  the  Severn  Estuary,  United  Kingdom  (data  reported  by  Odd  and  Rodger,  1986). 
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Lutoclines  have  been  found  to  persist  even  under  fairly  drastic  wave-induced  agitation,  as  a  result  of  the 
significant  negative  buoyancy  of  bottom  mud.  However,  under  highly  eroding  conditions,  may  decrease 
measurably  below  (Ross  and  Mehta,  1989)  and,  in  fact,  the  equality  between  and  is  valid  only  for 
conditions  under  which  the  sediment  settles  without  resusptension  or  upward  diffusion.  It  holds  reasonably  well  for 
weakly  eroding  conditions,  however. 

Fluid  mud  is  better  described  as  "fluid-like”  mud,  in  which  the  particles  are  largely  (but  not  always  solely) 
fluid-supported  (Smith  and  Kirby,  1989).  Thus,  in  general,  it  is  a  "quasi-suspension".  Since  the  bed  below  it  is 
essentially  particle-supported  the  density,  p„  is  the  value  of  p  below  which  inter-granular  contact  is  marginal.  Since 
if  left  at  rest  fluid  mud  will  dewater  to  form  a  bed,  its  existence  in  the  present  case  is  dependent  on  wave  agitation, 
which  prevents  a  reduction  in  the  pore  pressure,  and  thereby,  in  turn,  prevents  the  constituent  panicles  from 
developing  permanent  electro-chemical  bonds.  In  many  shallow,  wave-dominated  environments,  continued  wave 
action  causes  the  top  layer  of  mud  to  remain  fluidized.  For  instance,  in  Lake  Okeechobee,  Florida  wind  waves 
persistently  sustain  the  top  5-20  cm  thick,  organic-rich  mud  in  the  fluidized  state  (Kirby  et  al.,  1989;  Hwang,  1989). 

Since  the  fluid  mud/bed  boundary  is  typically  very  dynamic,  and  pore  pressures  difficult  to  measure  in  the 
field,  practical  definitions  have  been  used  to  determine  p,.  One  such  definition  is  based  on  the  approximate  empirical 
relationship  between  the  vane  shear  strength  of  the  soil,  and  p  or,  equivalently  the  solids  volume  fraction,  <i)-. 
r,=  ao(<f)-(f),)^“  ,  where  <t>,  is  solids  volume  fraction  corresponding  to  p„  and  Uq  and  are  sediment-specific 
coefficients  that  must  be  determined  experimentally.  Thus,  when  t,=0,  and  shearometric  evidence  has  been 
used  to  characterize  0,  as  the  critical,  space-filling  solids  volume  fraction  (James  et  al.,  1988).  For  the  mud  from 
Lake  Okeechobee,  Hwang  (1989)  obtained  (given  in  Pa)  ao=22.6,  1  and  (f>,=0.06;  the  latter  corresponding 

to  p,  =  1 ,065  kg/m^  Other  definitions  of  a  similar  nature,  but  ones  in  which  is  substituted  by  the  upper  Bingham 
yield  stress,  r,,  which  is  obtained  for  the  typically  pseudoplastic  stress-rate  of  strain  curve  for  mud  at  low  rates  of 
strain,  have  also  been  used  widely.  See  Mehta  (1991b)  for  a  brief  review  of  these  definitions.  In  recent  years  a 
revised  interpretation  of  as  the  critical  stress  at  which  plastic  yield  occurs  in  a  creep  test  has  been  proposed 
(James  et  al.,  1988;  Jiang,  1993).  This  test,  which  must  be  conducted  in  a  controlled-stress  rheometer,  provides 
a  direct  measure  of  t,,  thus  obviating  approximations  inherent  in  the  estimation  of  t,. 

The  most  commonly  considered  values  of  p„  and  p,  are  1,030  and  1,300  kg/m\  respectively,  although  other 
ranges  have  been  reported,  as  noted  in  Table  3  (amended  from  Ross  et  al.,  1987).  The  main  reason  for  the  observed 
variations  in  p^  and  p,  is  that  they  were  determined  under  different  hydrodynamic  conditions  and  for  different  muds. 

Cohesive  muds  having  concentrations  0  >  <^,  tend  to  exhibit  a  measurable  viscoelastic  response  to  wave 
forcing  of  stresses  below  the  critical  stress  for  plastic  yield.  The  viscosity,  p,  and  shear  modulus  of  elasticity,  or 
the  rigidity  modulus,  G,  of  typical  surficial  muds  found  in  the  shallow  marine  environment  range  from  O(IO')  to 
O(IO^)  Pa.s  and  Pa,  respectively.  These  ranges  indicate  that  such  muds  tend  to  be  highly  viscous  but  only  weakly 
to  moderately  elastic. 

By  way  of  an  elaborate  but  also  theoretically  more  accurate  definition  than  stated  above,  Foda  et  al.  (1993) 
consider  the  bed  to  be  an  elastic  medium,  and  "fluid"  mud  to  incorporate  a  transitional  viscoelastic  layer  between 
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Table  3.  Fluid  mud  density  and  corresponding  solids  volume  fraction  ranges 


Investigator(s) 

P' 

(kg/m^) 

Pt‘ 

(kg/m^) 

<t>,‘ 

Inglis  and  Allen  (1957) 

1,030 

0.018 

1,300 

0.182 

Krone  (196?>'’ 

1,010 

0.0061 

1,110 

0.067 

Wells  (1983) 

1,030 

0.018 

1,300 

0.182 

Nichols  (1985) 

1,003 

0.0018 

1,200 

0.121 

Kendrick  and  Derbyshire  (1985) 

1,120 

0.073 

1,250 

0.152 

Hwang  (1989) 

1,002 

0.047 

1,065 

0.060 

•Conversion  between  density  and  solids  volume  fraction  by  using  sediment  granular  density,  =  2,650 
kg/m^,  and  water  density  =  1,000  kg/m\  except  for  Hwang  (1989),  in  which  case  p,  =  2,140  kg/m\ 
•’Exclusively  based  on  laboratory  data.  The  other  definitions  rely  on  field  evidence,  although  they  are  not 
necessarily  based  solely  on  field  data. 

the  bed  and  a  viscous  layer  above.  Returning  however  to  the  more  simple  definition,  as  shown  in  Fig.  4,  typical 
shearometric  data  suggest  <i>,  to  be  in  the  range  of  0.05  to  0. 1 ,  below  which  G  is  comparatively  small  (James  et  al., 
1988).  <t>,  values  in  Table  3  on  the  other  hand  are  generally  higher,  with  the  exception  of  Krone’s  (0.067)  and 
Hwang’s  (0.060),  both  of  which  are  based  on  laboratory  data.  The  oti.>r  four  are  field-based,  and  technically 
somewhat  arbitrary.  At  the  same  time  however,  any  definition  of  d>,  based  on  G  (Fig.  4)  would  also  be  somewhat 
arbitrary,  unless  it  is  contingent  upon  a  quantitative  criterion  concerning  the  magnitude  of  G  for  defining  d>,. 

Since  the  dissipative  role  of  mud  is  largely  confined  to  the  top  compliant  layer  (Foda,  1989),  predictive 
hydrodynamic  models  that  assume  cohesive  mud  to  be  a  purely  viscous  fluid  seem  to  yield  simulations  of  surface 
wave  attenuation  of  acceptable  accuracy  for  field  application  (Jiang  and  Mehta,  1992).  In  other  simulations  an  elastic 


Figure  4.  Shear  rigidity  modulus  versus  solids  volume  fraction  for  sediment  SI  (after  James  et  al..  1988). 


component  has  been  added  (Jiang,  1993).  Less  dissipative  and  more  elastic  silty,  "fluid-like"  muds  have  been 
considered  to  be  poroelastic  (Yamamoto,  1983).  In  any  event,  it  is  essential  to  estimate  the  thickness  of  the  fluid 
(or  fluid-like)  mud  layer  as  a  precursor  to  hydrodynamic  modeling.  The  simplest  way  to  identify  layer  thickness 
is  in  terms  of  densities  and  p,.  However,  this  very  approximate  approach  does  not  explicitly  account  for  the 
dynamic  nature  of  the  fluid  mud  boundaries,  dependent  as  they  are  on  the  nature  of  wave  forcing.  Processes 
contributing  to  the  generation  of  fluid  mud  from  bed,  and  the  scour  of  the  fluid  mud-water  interface  leading  to  its 
entrainment  are  briefly  considered  next. 
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GENERATION  AND  EROSION  OF  FLUID  MUD 


3.1  FLUIDIZATION 

Under  progressive  waves  that  are  not  in  deep  water,  gradients  in  pore  pressure  can  cause  the  pore  fiuid  to  flow 
relative  to  the  soil  skeleton,  which  may  eventually  lead  to  rupturing  of  the  inter-particle  cohesive  bonds,  hence  a 
loss  of  effective  stress.  Starting  with  a  bed  at  rest,  the  increase  in  (wave-averaged)  pore  pressure,  u„,  with  time  at 
a  fixed  position  is  shown  schematically  in  Fig.  5,  in  which  o,,  =  hydrostatic  pressure  and  Au  =  wave-averaged 
excess  pore  pressure.  This  wave-averaged  effective  (normal)  stress,  o',  decreases  with  a  buildup  of  excess  pore 
pressure  until  soil  liquefaction  occurs  when  the  total  stress,  n=u,^=Au-fa|„  and  ff'=0  (see,  e.g.  Lamb  and 
Whitman,  1969).  The  wave-averaged  portrayal  of  pore  water  variation  is  an  evident  approximation  of  a  phenomenon 
that  actually  occurs  within  a  two-phased  medium,  in  which  properties  of  the  sedimentary  matrix  vary  spatially  at 
the  micro-fabric  scale,  even  when  the  medium  is  of  uniform  bulk  density  and  composition.  In  fact,  careful 
measurements  of  pore  pressure  oscillations  under  wave-forcing  in  silty  soils  suggest  that  transitional,  resonant 
amplification  of  the  pressure  amplitude  accompanies  an  episodic  increase  in  the  wave-mean  pore  pressure  at 
preferential  sites  (cavities)  within  the  soil  mass  (Foda  et  al.,  1991).  A  somewhat  similar  occurrence  is  reported  in 
a  clayey  soil  (Feng,  1992).  An  important  implication  of  this  observation  is  that  at  the  micro-fabric  scale,  soil 
liquefaction  is  initiated  in  a  spatially  and  temporally  inhomogeneous  manner. 

An  instantaneous  view  of  the  stress  profiles  in  the  mud-water  system  is  shown  schematically  in  Fig.  6.  When 
the  initial  bed  matrix  is  disturbed  by  wave  motion,  the  general  manner  in  which  fluidization,  when  considered 
synonymous  with  liquefaction,  occurs  depends  on  several  factors,  an  important  one  being  the  thickness  and  the 
degree  of  consolidation  of  the  mud  layer.  When  the  thickness  is  small  and  the  mud  is  unconsolidated  or  partially 
consolidated,  mud  displacement  tends  to  be  greater  near  the  rigid  bottom  below  compliant  mud,  than  near  the 
surface.  In  this  case  fluidization  can  conceivably  proceed  from  the  rigid  bottom  up.  In  a  comparatively  deep  mud 
layer,  or  a  consolidated  thinner  layer,  the  displacement  is  greatest  at  the  top  and  decreases  with  depth;  hence  in  this 


Figure  5.  Schematic  of  time-variation  of  pore  pressure  under  wave  action. 
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Figure  6.  Schematic  of  instantaneous  stress  profiles  in  a  water-mud  system. 


Figure  7.  Time-variation  of  effective  stress  at  three  elevations  (above  flume  bottom)  in  a  clay  bed  (kaolinite 
attapulgite)  subjected  to  1  Hz  progressive  waves  of  4  cm  amplitude  (after  Feng,  1992). 
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case  fluidization  tends  to  proceed  from  the  top  (Chou,  1989).  Formally  then  it  appears  that  the  manner  in  which 
fluidization  occurs  characteristically  depends  on  the  relative  water  depth,  kh,,  and  the  wave  Reynolds  number  in 

mud,  whj^/i/,,  where  k  =  wave  number,  w  =  wave  frequency,  h,  =  water  depth,  h,  =  mud  depth  and  v,  = 

kinematic  viscosity  of  mud.  For  small  values  of  kh,  and  tah/Zf,  boiiom-up  fluidization  should  be  favored.  For 
relatively  large  values  of  these  parameters,  ihe  reverse  process  should  be  facilitated,  provided  there  is  sufficient 
wave  energy  input  to  initiate  and  sustain  fluidization.  Imposing  a  non-uniform  bed  density  gradient,  with  density 
increasing  with  depth,  should  further  enhance  the  likelihood  of  fluidization  occurring  from  top,  since  in  this  case 
the  dynamic  pressure  gradient  tends  to  decrease  rather  markedly  with  depth  (Maa  and  Mehta,  1987). 

When  fluidization  proceeds  downward  from  the  mud/water  interface  an  equilibrium  thickness  of  the  fluid  mud 
layer  is  attained.  In  Fig.  7  the  fall  in  the  wave-averaged  effective  stress  with  time  at  three  selected  elevations  where 
time-series  of  total  and  pore  pressures  were  measured  below  the  interface  is  shown  from  a  preliminary'  test  (No. 
10;  see  Feng,  1992)  in  a  wave  flume.  The  16  cm  thick  bed  was  composed  of  an  aqueous  mixture  of  a  kaolinite  and 
an  attapulgite  in  equal  proportions  by  weight,  having  a  mean  density  of  1,170  kg/m^  The  water  depth  was  19  cm, 
and  the  1  Hz  forcing  wave  had  an  amplitude  of  4  cm.  The  bed  was  allowed  to  consolidate  fer  85  hours  before  test 
initiatioii.  Notice  the  increasing  time  required  for  fluidization  with  depth  below  the  interface. 

The  thickness  of  the  fluid  mud  layer,  d^t),  in  general  can  be  expressed  as 

d,  =  d,,-f(t)  (1) 

where  d^  is  the  equilibrium  value  of  df,  such  that  the  function  f(t)  -*  1,  as  t  -*  co.  Data  such  as  those  shown  in 
Fig.  7  were  used  to  determine  the  increase  in  d,  with  time  in  Fig.  8  from  three  tests  (Nos.  8,  9  and  10)  summarized 
in  Table  4  The  comparatively  short  durations  of  these  tests  were  insufficient  for  the  attainment  of  equilibrium 
depths,  df,.  In  any  event,  selecting  the  following  empirical  equation  for  fitting  the  data  points, 


Table  4.  Fluidization  depth  related  parameters 

Test  No. 

Wave 

Amplitude 

(cm) 

Consolidation  Period 
(hr) 

M, 

(cm) 

8 

2.0 

240 

18 

9 

2.8 

65 

50 

10 


4.0 


85 


."15 


15 


Figure  8.  Time-variation  of  the  depth  of  fluidization  in  three  tests  with  different  wave  amplitudes  and  consolidation 
periods.  Data  points  from  Feng  (1992), 


Figure  9.  Development  of  a  fluid  mud  layer  starting  with  a  "bed"  of  Hillsboro  Bay  mud  (Run  1;  see  Ross  and 
Mehta,  1990).  d,  is  the  vertical  distance  between  the  lutocline  and  the  a’  =  \  Pa  level.  Note  the  occurrence  of  an 
initial,  1.1  cm  thick  layer  of  fluid  mud  over  the  bed,  as  a  result  of  a  local  disturbance  of  the  bed  prior  to  test 
initiation.  The  observed  initial  increase  in  the  lutocline  elevation  is  attributed  to  advection  of  fluid  mud  to  the 
measurement  site. 
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we  note  that  the  curves  drawn  are  obtained  by  retaining  constant  values  of  the  coefficients  a|(=0.02l6)  and 
(=0.346)  for  the  time-dependent  function  {f(t)]  within  the  brackets,  while  varying  M,  (Table  4).  Thus  the 
influences  of  wave  amplitude  and  bed  consolidation  period  are  contained  wholly  in  M,.  Since  both  the  amplitude 
and  the  consolidation  period  were  varied  in  each  test  their  individual  effects  carmoi  be  gleaned;  however  we  note 
that  the  trend  of  increasing  rate  of  fluidization  corresponds  with  decreasing  bed  consolidation.  Thus  the  initial  state 
of  the  bed  seems  to  have  had  a  dominant  influence  on  the  rate  of  fluidiza,ion  in  these  panicular  tests.  Note  that 
according  to  the  form  of  Eq.  2,  M,=df5;  yet,  because  df,  was  not  actvially  attained,  M,  should  preferably  not  be 
set  equal  to  d^,  and  the  applicability  of  Eq.  2  should  be  restricted  to  the  experimental  durations  only. 

A  difficulty  with  tracking  wave-induced  fluidization  by  means  other  than  pore  pressures,  such  as  by  measuring 
accompanying  changes  in  bed  density,  is  that  conventional  methods  for  measuring  density,  e.g.  by  collecting  mud 
samples  for  subsequent  gravimetric  analysis,  are  often  too  coarse  for  the  degree  of  resolution  required  to  detect  the 
necessary  changes.  In  fact,  there  may  be  no  "significant"  change,  as  in  the  case  of  the  fluidization  test  shown  in 
Fig.  9,  in  which  the  time-variation  of  d^  is  shown  for  a  bed  composed  of  a  predominantly  montmorillonitic  estuarine 
mud  (from  Hillsboro  Bay,  Florida).  The  water  depth  was  31  cm  and  the  nominal  mud  thickness  12  cm.  The  bed, 
which  was  prepared  168  hr  before  test  initiation,  was  forced  by  a  1  Hz  wave  of  3  cm  amplitude.  The  boundary 
between  the  fluid  mud  layer  and  the  cohesive  bed  is  defined  by  the  locus  of  the  elevation  correspondin'  to  a'  =  1 
Pa,  a  very  small  value.  During  the  first  hour,  vertical  profiles  of  mud  density  were  measured  (by  gravimetric 
analysis  of  samples  withdrawn  from  the  side  of  the  experimental  flume  via  tubes),  as  shown  in  Fig.  10.  Arrows  1 
through  4  mark  the  positions  of  the  a'  =  I  Pa  boundary  in  Fig.  9.  Note  that  between  levels  1  and  4  no  measurable 
change  in  density  occurred  (Ross  and  Mehta,  1990). 

Sensitive  measures  of  "soil  softening"  by  waves  include  such  lumped  parameters  as  the  undrained  shear 
strength  and  the  shear  modulus  or  rigidity  modulus,  G,  of  the  bed  (Thiers  and  Seed,  1968;  Schuckman  and 
Yamamoto,  1982).  V/i;hin  'he  bed  of  kaolinite  and  attapulgite  mentioned  earlier,  and  having  a  pre-test  consolidation 
period  of  20  hr,  a  specially  designed,  miniature  shearometer  was  installed  to  record  the  change  in  G  with  time  under 
continued  wave  action  (Williams  and  Williams,  1992).  The  device  measured  the  speed  of  propagation,  V,  of  a 
1.8  kHz  shear  wave  generated  and  detected  by  transducers  incorporating  piezo-ceramic  "bimorph"  elements  bonded 
to  miniature  steel  plates  which  served  as  the  generating/detecting  surfaces  (length  12  mm;  width  4  mm;  thickness 
100  /tm).  At  this  high  shear  wave  frequency,  V  was  practically  instantaneous  in  terms  of  its  variation  over  the 
experimental  time  scale.  Given  the  initial  speed  V(0),  the  corresponding  ratio,  V(t)/V(0),  is  plotted  against  time  in 
Fig.  11.  Notice  the  "exponential"  fall  in  the  speed,  which  is  commensurate  with  a  corresponding  increase  in  df  such 
as  that  observed  in  Fig.  8.  The  best  fit  equation  is 

=  e  (3) 

V(0) 

where  aj =0.0216  ana  =0.346.  Recognizing  that  the  forms  of  the  time-dependent  terms  in  Eqs.  2  and  3  are 
analogous,  and  that,  perhaps  fortuitously,  «,=«,  and  /3,=/32,  eliminating  time  between  the  two  equations  leads  to 
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Figure  10.  Mud  density  profiles  from  Run  1  (Fig.  9).  Arrows  indicate  bed  level  (a'  =  1  Pa)  at  different  times  at 
test  initiation  (level  1),  and  following  test  initiation  (after  Ross  and  Mehta,  1990). 


Figure  1 1 .  Time-variation  of  relative  shear  wave  velocity  in  a  bed  of  attapulgite  -t-  kaolinite  of  20  hr  consolidation 
subjected  to  a  1  Hz  progressive  water  wave  of  2  cm  forcing  amplitude  (after  Williams  and  Williams,  1992). 
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d,  =  M,  1 


vcoV’ 

V{0)_ 


(4) 


which  indicates  a  seemingly  unique  dependence  of  df  on  V  for  a  given  bed.  For  a  given  shear  wave  forcing 
frequency,  G  depends  on  V  as  well  as  on  the  phase  angle  between  the  shear  stress  amplitude  and  the  resulting 
amplitude  of  strain  (Williams  and  Williams,  1992).  Thus  the  coefficient  M,  essentially  contains  the  influence  of  the 
phase  angle  in  an  unquantifiable  way.  Nevertheless,  Eq.  4  demonstrates  the  dependence  of  d(  on  a  soil  dynamical 
parameter  that  intrinsically  reflects  the  physical  state  of  mud. 


3.2  ESTIMATION  OF  FLUIDIZATION  THICKNESS 

Research  on  the  estimation  of  the  fluidization  thickness  apparently  has  largely  focussed  on  the  equilibrium 
value,  dfj,  as  opposed  to  d^t).  Thus,  for  example,  the  linkage  between  mud  rheology  and  df^  was  explored  by  Chou 
(1989),  by  considering  bottom  mud  to  be  generally  viscoelastic.  Onset  of  fluidization  was  defined  in  terms  of 
critical  values  of  shear  strain  related  to  elastic  and  viscous  responses  of  the  bed.  This  criterion  is  mathematically 
different  from  that  for  liquefaction  (i.e.  a’  =0).  Yet,  the  two  characterize  the  same  phenomenon.  As  noted  in  Section 

4.4.3,  the  constitutive  equation  for  a  viscoelastic  material  is,  T,,  =G*Ej,,  where  T,,  =  devialoric  component  of 
stress,  Eji  =  deviatoric  component  of  strain,  subscripts  i,j  denote  directions  and  G*  is  the  complex  shear  modulus. 

By  definition,  G*=G'-iG",  where  G'  =  storage  modulus  (associated  with  the  elastic  component)  and  G"  =  loss 
modulus  (associated  with  the  viscous  component).  For  a  Voigt  material,  G'=G  and  G"=/iw,  where  n  =  viscosity, 
G  =  shear  modulus  of  elasticity  and  u  =  forcing  frequency  (Barnes  et  al  .,  1989). 

The  behaviors  of  G'  and  G"  as  functions  of  the  amplitude  of  applied  oscillatory  strain,  y\  were  examined 
by  Chou  (1989)  in  a  controlled-strain  rheometer.  Figure  12  shows  the  results  for  a  kaolin.  Note  the  importance  of 
elasticity  at  low  strains  when  viscous  loss  was  practically  nil.  With  increasing  strain  the  clay  behaved  increasingly 
as  a  viscous  material.  Results  such  as  these  were  used  to  develop  the  nomograms  of  Figs.  13a, b  for  G’  and  G”  as 
functions  of  strain  amplitude,  7°  and  relative  mud  density  (p/p»).  The  forcing  frequency  was  held  constant  at  1.5 
rad/s.  However,  in  experimenting  with  changing  the  forcing  frequency  it  was  found  that  the  influence  of  7  was 
much  more  pronounced  than  that  of  w  in  controlling  G'  and  G”.  Depending  on  the  value  of  7%  the  domain  is  seen 
to  be  subdivided  into  three  regimes  -  viscous,  viscoelastic  and  elastic;  7„  and  are  values  of  the  strain,  7,  that 
define  the  appropriate  boundaries. 

The  G',  G"  nomograms  enable  an  empirical  determination  of  the  thicknesses  of  the  viscous  and  viscoelastic 
layers  in  an  otherwise  elastic  bed.  Chou  (1989)  considered  the  water-mud  system  shown  schematically  in  Fig.  14, 
and  calculated  the  equilibrium  thicknesses  of  the  fluid-like  viscous  layer  and  the  soft  viscoelastic  layer  as  functions 
of  the  wave  height  H  =2a,  where  a  is  the  forcing  wave  amplitude.  The  procedure  involved  the  use  of  a  wave-mud 
interaction  model  in  conjunction  with  the  G',  G”  nomograms.  For  a  Voigt  viscoelastic  mud  the  model  becomes 
analogous  to  that  of  Maa(1986).  The  model  was  used  iteratively  to  calculate  strain  amplitudes  at  different  depths 
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Figure  12.  Storage  modulus,  G'(*),  and  loss  modulus  G"(o),  versus  strain  amplitude,  7  %  for  K2  kaolin  at  a  salinity 
of  35  ppt  and  a  solids  concentration  of  4 1  % .  w  =  1 .5  rad/s.  Lines  are  drawn  to  enhance  data  trends  (after  Chou, 
1989). 
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Figure  13.  Contours  of  storage  modulus,  G'  in  (a)  and  loss  modulus,  G",  in  (b)  for  K1  kaolin  at  a  salinity  of  35 
ppt  (after  Chou,  1989). 
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Figure  14.  A  four-layer  viscoelastic  model  (adapted  from  Chou.  1989). 


WAVE  HEIGHT  (m) 


Figure  15.  Thickness  of  fluid  mud  (for  a  clayey  sediment)  as  a  function  of  wave  height.  H.  assuming  a  water  depth 
of  5  m,  a  wave  frequency  of  0,1  Hz  (adapted  from  Chou.  1989). 
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as  iimctions  of  the  wave  amplitude,  and  hence  the  depths  h,  and  hj.  An  example  of  the  fluidized  or  fluid-like  layer 
thickness,  hj-h,,  as  a  ftmction  of  wave  height  is  shown  in  Fig.  15. 


3.3  SURFACE  EROSION 

Since  the  action  of  waves  over  a  cohesive  bed  leads  to  a  fairly  rapid  development  of  the  fluidized  mud  layer, 
and  since  fluidization  can  be  considered  to  be  a  mode  of  entrainment,  it  is  sometimes  reasonable  to  consider  the 
interface  between  fluid  mud  and  the  bed  to  be  the  appropriate  bed  "surface”.  This  assumption  however  becomes 
increasingly  tenuous  with  increasing  thickness  of  the  fluid  mud  layer.  Accurate  models  for  sediment  transport  must 
therefore  consider  bed  fluidization  as  a  process  that  is  distinct  from  mud  surface  erosion.  Nevertheless,  in  a  large 
number  of  practical  applications  the  use  of  stress-based  equations  for  the  rate  of  bottom  erosion,  akin  to  those 
obtained  for  uni-directional  flow  situations,  seems  to  yield  acceptable  answers  with  regard  to  the  entrainment  of 
particles  into  the  water  colunm. 

Expressions  for  the  rate  of  erosion  have  been  derived  almost  exclusively  from  laboratory  studies.  Exceptions 
include  "back  calculated"  rates  using  field  data  and  numerical  models  for  suspended  sediment  transport  (e  g.  Sheng 
et  al.,  1986).  Some  laboratory  erosion  rate  expressions  are  listed  in  Table  5.  The  expression  of  Alishahi  and  Krone 
(1964)  is  based  on  data  on  the  time-variation  of  the  depth-mean  suspension  concentration  during  erosion  reported 
by  the  investigators.  Experimental  data,  in  general,  typically  conform  to  the  relationship 


e  = 


(5) 


where  e  =  rate  of  erosion,  =  peak  value  of  the  cyclic  bed  shear  stress,  Mj  =  erosion  rate  constant,  =  erosion 
resistance  and  a,  =  empirical  constant.  Values  of  a,  in  Table  5  range  from  0.95  to  1.82.  When  0:3=1,  Eq.  5 
becomes  analogous  to  the  corresponding  expression  for  cohesive  bed  erosion  under  steady  or  quasi-steady  flows  (see 
e.g.  Mehta,  1988): 


e  =  M, 


(6) 


where  M3  is  the  rate  constant  and  r,  is  defined  as  the  erosion  shear  strength.  In  the  laboratory  setting,  r,  can  be 
determined  as  a  function  of  depth  below  the  bed  surface  via  a  procedure  involving  layer-by-layer  erosion  of  the  bed 
(Parchure  and  Mehta,  1985).  Maa  (1986)  used  a  similar  procedure  to  determine  Tr,  and  thus,  in  a  sense,  established 
a  correspondence  between  r,  and  tr.  Furthermore,  the  range  of  0.15  to  0.39  Pa  reported  for  tr  in  Table  5  is 
generally  of  the  same  order  of  magnitude  as  that  for  t,  (see  e.g.  Mehta,  1988).  However,  it  must  also  be  recognized 
that  the  time-change  in  the  properties  of  the  bottom  mud  due  to  wave  action  tends  to  influence  Tr  in  a  significant 
way  (Maa,  1986).  Thus  the  value  of  Tr  used  in  solving  Eq.  5  must  be  inherently  a  characteristic  time-mean  value, 
unless  the  variation  of  Tr  with  time  is  explicitly  included.  However,  given  the  usual  paucity  of  data  on  Tr  and 
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empirical  constant;  U  =  aw/sinh(kh,)  is  the  amplitude  of  the  bottom  oscillatory  velocity,  where  k  =  wave  number  and  h,  =  water  depth;  5^  =  (i'/2w)''’  is 
the  wave  boundary  layer  thickness  where  v  =  kinematic  viscosity  of  water;  f(z)  =  bed  shear  strength  as  a  function  of  depth  (z);  Test  2  (148  hr)  had  a  longer 
consolidation  period  than  Test  1  (38  hr). 


unquantifiable  effects  of  bed  changes  on  particulate  entrainment,  the  use  of  time-varying  is  generally  not 
warranted.  It  is  also  noteworthy  that,  depending  on  the  way  in  which  waves  influence  the  bed,  Tr  may  be  lower  or 
higher  than  r,  (Mimura,  1993).  In  Fig.  16,  Tr  is  observed  to  be  lower  than  t,  for  a  bed  of  kaolinite. 

In  Fig.  17,  a  typical  relationship  in  the  form  of  Eq.  5  based  on  flume  measurements  is  shown  (Maa  and  Mehta, 
1987).  Data  smearing  of  the  degree  observed  is  typical  of  such  plots,  and  is  at  least  in  part  due  to  the  uncenainty 
in  estimating  Tr.  Furthermore,  it  is  rather  self-evident  that  plots  such  as  these  mask  influences  that  may  depend  on 
multiple  causes  and  associated  feedbacks.  For  example,  Jackson  (1973),  who  conducted  laboratory  tests  on  wave 
resuspension  of  a  mixture  of  silt,  clay  and  some  fine  sand  on  a  sloping  beach,  observed  rapidly  increasing  suspended 
sediment  concentration  in  the  first  hour,  reaching  a  maximum  after  the  second  hour  and  decreasing  subsequently. 
A  possible  explanation  for  the  decrease  in  concentration  is  that  the  bed  became  less  rigid  with  time  under  wave 
action,  which  in  turn  resulted  in  greater  energy  dissipation  and  reduced  the  shear  stress  at  the  interface.  Jackson 
also  noted  that  for  relatively  small  orbital  velocities  (12-15  cm/s),  there  was  very  little  erosion,  whereas  for 
velocities  greater  than  about  20  cm/s  the  surface  erosion  process  was  replaced  by  mass  erosion,  with  comparatively 
much  larger  rates  of  sediment  entrainment,  a  behavior  also  observed  under  steady  flows  (Mehta,  1991c). 

3.4  BED  RECOVERY 

Since  the  experimental  data  presented  in  Section  3.2  suggest  that  the  loss  of  effective  stress  is  primarily  due 
to  the  rupmring  of  inter-particle  bonds  without  a  significant  change  in  the  density  or  void  ratio,  it  follows  that  bed 
recovery  should  be  comparatively  rapid  in  the  laboratory  setting  after  wave  action  ceases.  In  test  No.  9  of  Feng 
(1992)  mentioned  earlier,  wave  action  was  maintained  for  7.5  hr,  following  which  effective  stresses  at  different 
elevations  above  the  flume  bottom  were  determined  from  pore  pressure  profiles  4  hr  later.  Observe  in  Table  6  that 
the  percent  recovery,  as  defined,  varied  from  0  to  77,  with  a  mean  of  36.  The  non-uniformity  in  percentages 
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Figure  16.  Influence  of  waves  on  shear  resistance  to  erosion  of  kaolinite  beds  in  flumes  (after  Mehta,  1989). 
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Figure  17.  Normalized  rate  of  erosion  versus  normalized  excess  shear  stress  for  a  bay  mud  (after  Maa  and  Mehta, 
1987). 

Table  6.  Bed  recovery  in  Test  No.  9  of  Feng  (1992) 


Elevation 

(cm) 

t=0 

a'(Pa) 

t=7.5  hr 

t  =  11.5hr 

Recovery 

12.0 

39 

8 

8 

0 

9.5 

53 

0 

41 

77 

7.5 

67 

7 

35 

47 

5.1 

83 

53 

59 
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presumably  reflects  local  variability  in  the  mud  properties  at  the  sites  of  pressure  measurements.  In  general, 
however,  with  the  exception  of  the  data  point  at  12  cm  (4  cm  below  the  interface),  recovery  is  observed  to  slow 
down  with  increasing  depth.  Overall,  since  the  attapulgite  +  kaolinite  mixture  was  thixotropic,  the  rate  of  recovery 
is  likely  to  be  related  in  some  way  to  the  rate  of  gelling.  Thus,  for  example.  Day  and  Ripple  (1966)  reported  a 
recovery  period  of  14  hr  for  an  initially  sheared  K-montmorillonite  based  on  the  measurement  of  suction  in  a 
tensiometer.  Typically,  gelling  is  complete  in  about  a  day. 
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ENERGY  DISSIPATION  AND  ELEMENTS  OF  MODELING 


4.1  WAVE  ENERGY  DISSIPATION  OVER  NON-RIGID  BEDS 

The  conservation  equation  for  wave  energy  within  the  water  column  is: 

^  +  V,<qE)^e„  =0  (7) 

where  E  =  mean  wave  energy  per  unit  surface  area,  =  horizontal  gradient  operator,  =  wave  group  velocity 
and  Cd  =  rate  of  dissipation  of  wave-mean  energy  per  unit  area.  Constitutive  models  proposed  for  the  rate  of  energy 
dissipation  within  mud  generally  fall  into  two  groups:  1)  models  which  are  specific  to  the  assumed  description  of 
mud  rheology,  and,  2)  models  which  are  non-specific  with  regard  to  rheology.  Some  of  the  proposed  expressions 
for  fo  aic  listed  in  Table  7.  In  the  first  group,  cohesive  mud  has  been  considered  either  as  a  viscous  fluid,  or  as 
a  viscoelastic  or  poroelastic  medium.  Some  investigators,  e.g.  Shibayama  et  al.  (1990),  have  also  considered 
cohesive  mud  to  be  a  Bingham  viscoplastic  in  which,  at  stresses  less  than  the  yield  value,  the  material  is  treated  as 
a  viscoelastic,  whereas  at  higher  stresses  it  is  a  fluid.  Silty  muds,  considered  to  be  poroelastic,  incorporate  energy 
loss  by  Coulomb  friction  between  clay  particles  (e.g.  Yamamoto  and  Takahashi,  1985). 

In  the  second  group,  the  expression  of  Tubman  and  Suhayda  (1976)  requires  coefficients  M  and  0  which 
depend  on  the  character  of  the  water  pressure  wave  near  the  mud  surface  and  that  of  the  induced  mud  surface  wave. 
Since  these  waves  depend  on  the  solution  of  the  hydrodynamic  problem  with  an  assumed  constitutive  relationship 
for  mud,  M  and  <t>  must  be  measured  directly,  or  obtained  indirectly  via  model  calibration  using  experimental  data. 
Likewise,  the  expression  of  Schreuder  et  al.  (1986)  involves  the  mean  shear  strain  rate  as  an  explicit  input 
parameter.  Appendix  A  provides  additional  information  on  the  bases  of  the  eo  expressions. 

4.2  A  MODELING  PERSPECTIVE 

The  response  elements  shown  in  Fig.  1,  namely  wave  attenuation,  mud  motion  and  surface  erosion,  do  not 
easily  lend  themselves  to  a  unified  theoretical  treatment,  since  the  first  two  can  be  examined,  at  least  in  the 
simplified  treatment,  by  assuming  bottom  mud  to  be  a  constant  density,  incompressible  continuum,  whereas 
particulate  entraiiuncnt  at  the  mud-water  interface  requires  separate  consideration  in  terms  of  the  change  in  the 
density  of  the  overlying  water  layer  due  to  a  fluid  stress-dependent  vertical  mass  flux  at  the  interface.  Even  under 
the  most  severe  wave-induced  agitation  of  bottom  mud,  the  lack  of  a  competent  diffusive  mechanism  to  transport 
the  sediment  mass  from  the  bottom  to  the  top  of  the  water  column  usually  precludes  the  development  of  high 
suspension  concentrations  in  the  surficial  waters;  most  of  the  material  remains  close  to  the  bed,  approximately  within 
the  wave  boundary  layer.  In  turn,  the  depth  of  scour  of  bottom  mud  is  typically  sufficiently  small  to  make  it 
practicable,  at  least  to  first  order  accuracy,  to  assume  water  and  mud  thicknesses  to  be  time-independent.  Under 
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Table  7.  Some  energy  dissipation  rate  expressions  due  to  non-rigid  bed  mechanisms 


Mud  property* 


Group  1 

Mud  as  viscous 
fluid’ 


Mud  as  viscous 
fluid^ 

Mud  as  viscous 
fluid^ 


Expression” 


m 

1 

2*'2 

e’'‘>V-gk,)’ao^ 

4v5' 

0) 

^D1  +  £d2  ^D3 

^04  +  £dS 


Investigator(s) 


Dean  and  Dalrymple  (1984) 


Dalrymple  and  Liu  (1978) 


Feng  (1992) 


Mud  as  Voigt 
body^ 


Coulomb  friction 
in  bed* 


Coulomb  friction 
in  bed* 


r  1 

Pig' 

4 

Pig' 

4 

<j  A,  Im 


<w^(gkcoshlch,  -w^siiihkh,) 
gk(gksinhkh,  -w^coshkh,) 


MacPherson  (1980) 


Yamamoto  and  Takahashi 
(1985) 


Yamamoto  and  Takahashi 
(1985) 


Group  2 
Mud  wave’ 


Mud  as 

viscoelastic 

medium* 


Fi8 

2 


wMsin<i)a’ 

cosh’kjh. 


(ft'  <7>’hj) 


Tubman  and  Suhayda  (1976) 


Schreuder  et  al.  (1986) 


•See  Appendix  A  for  explanatory  notes  1  through  8. 

”Notation:  1  =  water  layer;  2  =  mud  layer;  p  =  density;  v  =  kinematic  viscosity;  w  =  wave  frequency;  k  = 
complex  wave  number  =  k^-i-iki,  where  Iq  =  wave  attenuation  coefficient,  k,  =  wave  number  and  i=(-l)''’;  h  = 
depth;  g  =  acceleration  due  to  gravity;  a  =  local  wave  amplitude;  Sq  =  refeience  wave  amplitude;  ep,  through 
=  rate  of  energy  dissipation  terms;  A,  =  wave  amplitude  ratio;  Im  =  imaginaiy  part  of  a  complex  variable;  5  = 
Coulomb  specific  loss;  G  =  dynamic  shear  modulus;  ko  =  wave  number  for  a  rigid  bottom;  osq  =  2Gko’/p2;  M 
=  proportionality  constant  between  the  amplitudes  of  the  mud  wave  and  the  wave-induced  bottom  pressure  wave 
;  <^=180®-®,  where  9  =  phase  angle  between  the  crest  of  the  bottom  pressure  wave  and  the  crest  of  the  mud  wave; 
fjt'  =  dynamic  viscosity;  <7>  =  average  value  of  shear  strain  rate. 


24 


these  assumptions,  solutions  developed  for  problems  involving  stratified,  non-Newtonian  fluid  mechanics  have  been 
applied  extensively  to  model  wave-mud  interaction.  Some  of  the  previous  developments  are  summarized  in 
Appendix  B. 

In  what  follows,  the  formulation  of  the  linearized  hydrodynamic  problem  for  wave  attenuation  and  mud  motion 
is  considered  in  conjunction  with  Newtonian  viscous  and  linear  viscoelastic  descriptions  of  mud  rheology. 
Progressive,  non-breaking  waves  are  treated  for  their  propagation  within  shallow  water  and  non-shallow 
(intermediate  depth)  water  conditions.  Following  this  development,  simple  bases  for  modeling  erosion  are  discussed 
briefly.  These  presentations  are  for  informative  purposes  only;  for  details  the  referenced  works  in  Appendix  B  must 
be  consulted. 


4.3  BASIC  CONSERVATION  EQUATIONS 

For  continuous  media,  the  governing  equation  of  continuity  is: 


=  0 

Dt  '^ax, 


(8) 


where  pfXpt)  =  density,  U|(Xi,t)  =  velocity  and  subscript  i  denotes  direction.  For  an  incompressible  fluid  Dp/Dt=0, 
and  the  continuity  equation  is  simplified  to 


dx. 


(9) 


The  corresponding  equation  of  momentum  is: 


p.P!!i=pb 

^  Dt  ^  '  dx. 


(10) 


where  bj(Xj)  =  body  force,  Tj|(X|,t)  =  external  surface  stresses,  and  subscripts  i  and  j  denote  directions  in 
conjunction  with  subscript  ji,  which  implies  a  second  order  tensor.  Thus,  j  =  direction  normal  to  the  surface  and 
i  =  stress  direction.  Because  gravity  is  the  only  body  force  of  concern  here  and  as  we  set  the  z-tiirection  vertically 

upward,  then  b|=b2=0  and  b3=-g.  The  stress,  T„,  can  be  expressed  as  the  sum  of  the  deviatoric  component, T„ 
and  the  spherical  component,  .iTu,6j|  =  -paj„  where  T^  is  the  sum  of  the  diagonal  terms  in  the  stress  tensor,  p  is 
the  mean  normal  pressure  and  =  kronecker  delta;  for  j=i,  5j|=l,  and  for  j?ii,  6j,=0.  Thus, 


DUi  dp  dT„ 

'’-dT  ■  ■  ar  ar 


(11) 
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4.4  CONSTITUTIVE  EQUATIONS 

4.4.1  Fluid  (incompressible) 

The  constitutive  equation  for  the  fluid  is; 


Tj,  =  2/iiDj, 


where  p  =  viscosity,  and  the  rate  of  strain,  Dj,  =0.5[(du,/dx,)+(du/dXj)].  Thus,  the  momentum  equation  becomes 


Du,  Sp  d\ 


dx,  dx.dx. 


4.4.2  Elastic  Material 

The  constitutive  equation  for  elasticity  is: 

T,;  =  2GE,i  (14) 

where  G  =  modulus  of  elasticity,  and  the  strain  Eji  =0.5[(dL,/3X|)+(dLs/dXj)],  where  L„  Lj  are  displacements.  Thus 
the  momentum  equation  becomes 

=  pSjb  -  iE  +  G-!!h.  (15) 

Dt  Dt^  ’’  ■  dx.  dx^dx^ 


4.4.3  Viscoelastic  Material 

The  general  constitutive  equation  for  a  linear  viscoelastic  medium  is 


M  N 

_  O  .  1- 


rM)  Ol  5.K)  Ol 


where  subscripts  r  and  s  denote  orders  of  partial  differentiation,  p,  and  q,  are  coefficients  related  to  the  viscoelastic 
model,  and  M  and  N  are  the  specified  maximum  differential  orders  of  the  chosen  model.  Under  cyclic  loading,  we 


consider 


Tji  =  Toexpf-i(wt)I 


Ej,  =  Eoexp[-i(«t-5)] 


26 


Note  that  for  a  viscoelastic  medium,  from  Eq.  19  we  can  also  define  T,,  =G*E|,.  where  G*=G'-iG”,  where  G' 
IS  the  elastic  energy  storage  modulus  and  G”  is  the  viscous  energy  dissipation  modulus.  Thus  it  can  be  shown  that 
for  the  Voigt  model,  for  example,  G'  -G  and  (Barnes  et  al..  1989). 

For  a  viscoelastic  material  -  Maxwell  model: 


2G 


;  = 


2Gw 


(26) 


^2 


Next,  as  illustrations  of  more  complex  models,  for  the  three -element  viscoelastic  models,  constructed  from  Voigt 
and  Maxwell  elements,  shown  in  Figs.  18a, b; 

Model  (a): 


2G, 

2i 

2G, 

•> 

f 

\  2 

;  P”  = - 

\  2 

G,  +G, 


(27) 


This  model  has  been  used  in  describing  the  consolidation  behavior  of  soils  (Keedwell,  1984).  Jiang  (1993)  used  it 
in  his  wave-mud  interaction  model  described  in  Section  4.5. 

Model  (b): 


2Gj 

2l 

2 

w 

A 

G, 

2  *  ^ 

G, 

2 

Ih 

(28) 


For  any  other  (linear)  viscoelastic  model,  the  complex  dynamic  viscosity,  n*=n'  +  can  be  obtained  through 
Eq.  20. 


4.5  MODELS  FOR  WATER  WAVES  OVER  MUD  (X-Z  PLANE) 

4.5.1  Inviscid,  Shallow  Water  over  Mud 

A  general  definition  diagram  for  the  wave-mud  system  is  shown  in  Fig.  19.  In  the  case  considered  first,  watei 
is  assumed  to  be  inviscid  and  mud  to  be  viscoelastic.  Under  the  shallow  water  assumption,  the  following  analytical 
solution  of  the  linearized  problem  is  obtained  based  on  the  works  of  Gade  (1957,  1958)  and  Jiang  (199?). 

For  water: 

The  continuity  Eq.  9  is  integrated  from  the  water-mud  interface  to  the  water  surface; 


r, 

f 

du, 

dw, 

L 

dx 

dz 

dz=0 


(29) 
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Figure  18.  Three-element  viscoelastic  models:  a)  Voigt  element  (represented  by  2G2)  modulated  by  additional 
elasticity;  b)  Maxwell  element  (represented  by  2^2.  2G2)  constrained  by  elastic  modulus  represented  by  2G,. 
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Figure  19.  Water-mud  system  in  the  x-z  plane. 


in  which  water  velocity,  Uj,  is  invariant  over  water  depth,  and  at  the  surface  and  the  interface: 


w 


ar, 

dt 


W 


N-fi 


d\ 


(30) 


Therefore,  the  continuity  equation  for  water  becomes: 


'  dx  dt 


(31) 
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FromEq.  13,  omitting  the  viscous  terms,  forinviscid  "water",  using  the  hydrostatic  assumption 

for  pressure,  p=Pig({’j-z),  in  shallow  water  and  omitting  the  non-linear  teniK,  u^du/SXj,  the  momentum  equation 
for  water  (in  the  x-direction)  becomes 


au, 

'  at 


ax 


^■®ar 


(32) 


For  mud: 

Equation  9  is  integrated  from  the  mud  bottom  to  the  water-mud  interface  using: 


w 


-h,-h, 


=  0 


w 


at 


Then  the  continuity  equation  for  mud  becomes 


-N-f. 


auj 

ax 


dz 


.S.D 

at 


(33) 


(34) 


From  the  momentum  Eq.  22,  using  the  hydrostatic  pressure  assumption  and  omitting  the  non-linear  terms,  we 
obtain: 


(P2-Pi)g 


ax 


(35) 


where  pj  =  mud  density.  The  horizontal  diffusion  term,  p'fa^j/ax^),  has  been  ignored  in  comparison  with  the  last 
term  on  the  right  hand  side  (in  Eq.  35)  corresponding  to  vertical  diffusion.  Eq.  35,  in  general,  is  applicable  to  small 
strain  oscillatory  forcing. 

The  following  boundary  conditions  are  next  introduced; 

At  the  free  surface,  z=0: 

f,(0,t)  =  Acos(«t)  (3<5) 


where  A  =  wave  amplitude. 

At  water-mud  interface,  z=-h,: 

au^(x,-h,.t)  ^  ^  (37^ 

az 


At  the  mud  bottom,  z=-(h,-hh2): 

Uj(x,-(h,+h,).t)=0  (38) 

Under  cyclic  wave  loading,  the  following  harmonic  assumptions  occur. 
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fi  =  f,  exp[i(kx-wt)]  (39) 

fj  =  f j  exp[i(kx  -tut))  (40) 

u,  =u,  exp[i(kx-tut)]  (41) 


Uj=Uj(z)  exp[i(lcx-<ut)] 


(42) 


where  the  amplitudes  denoted  by  are  independent  of  x,  z  and  t  except  in  Eq.  42,  wherein  Uj  varies  with  z.  Next, 
substituting  the  boundary  conditions  into  the  governing  equations  (31,  32,  34,  35)  and  after  some  manipulation  the 
following  solution  for  the  complex  wave  number,  k,  is  obtained: 


1+^r- 


h, 


N2 


-4r^r 


hi 

2r-3r 

“i 


(43) 


where  r=(p2-/>i)/P2:  Froude number,  F,=<u(h,/g)‘'^;  r=l-[tanh(mh2/h,)/(mhj/h|)];  m=(Rt  and  Reynolds  number, 

R, = wh,Vi/  *  (v  *  *  /pj) .  Note  that  k = k, + ikj,  where  k,  is  the  wave  number  and  Iq  is  the  wave  attenuation  coefficient . 

Thus  the  wave  amplitude  at  a  distance  x,  a,=a(^xp(-k|X),  where  ao  is  the  amplitude  at  x=0.  Solutions  for  the  other 
variables  are: 

=  Aexp[i(kx-<ut)]  (44) 

=  A 


kh, 


exp[i(kx-cut)) 


(45) 


kh, 

u,  =  wA — I  exp(i(kx-tut)] 
F' 


(46) 


kh, 

Uj  =  cuA  ' 


1-cosh 

z+(h  +hj) 
m  ' 

+tanh 

m— 

sinh 

z+(h|+lt,) 

m 

h, 

h,] 

h, 

(47) 


1-r 


kh, 

'W 


exp{i(kx-cjt)J 


Finally,  the  shear  stress  in  the  water  layer,  hence  at  the  interface,  is  zero  because  of  the  inviscid  assumption, 
and  in  mud  layer  the  shear  stress  is  obtained  from: 
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T(X.Z,t)  =  fl' 


-.^sinh 

2+(h,+hj) 

m 

h, 

+  — tanh 
hi 

"'tt 

h, 

cosh 

z+(h  +h,) 

r  ^  IJ 

kh, 

f  kh,]^ 
F,' 

expli(kx-<rt)] 

(48) 


Without  invoking  the  shallow  water  assumption,  Hsiao  and  Shemdin  (1980)  and  MacPherson  (1980) 
independently  investigated  the  problem  of  inviscid  water  waves  propagating  over  a  linear  viscoelastic  mud,  and 
obtained  analytical  solutions  in  which  the  wave  number  is  implicitly  expressed.  Suhayda  (1986)  extended 
MacPherson’s  results  to  a  non-linear  viscoelastic  bed,  while  Mei  and  Liu  (1987)  considered  a  viscoplastic  bed.  Foda 
(1989)  examined  a  stratified  elastic  bed  and  Feng  (1992)  and  Isobe  et  al.  (1992)  set  up  their  models  considering  mud 
to  be  a  non-Newtonian  (power-law)  fluid.  However,  as  noted,  since  in  these  models  the  viscosity  of  water  was  not 
included,  the  shear  stress  on  water-mud  interface,  which  governs  particulate  entrainment  at  the  interface,  was 
neglected. 


4.5.2  Viscid  Fluid  (Water)  over  Viscid/Viscoelastic  Mud 

Considering  the  viscosity  of  water  and  non-shallow  water  conditions,  the  governing  equations  can  be  set  up 
as  follows; 

From  Eq.  9  the  continuity  equation  is: 


da  dW;  ^ 
_ !  +  _ i=0 


(49) 


dx  dz 

where  u  and  w  are  velocities  in  the  x-  and  z-directions,  respectively,  and  i  =  1,2  denote  the  upper  (water)  and  lower 
(mud)  layers. 

From  Eq.  22  for  a  viscid  fluid,  omitting  the  non-linear  terms,  the  momentum  equation  coWponents  are; 


du,  1  dPj  d^i,  d^i- 

dt  Pi  dx  '  dx^  dz^ 


(50) 


dw,  1  dP,  d^v,  d^w, 
dt  p.  dz  '  dx^  dz^ 


(51) 


where  the  dynamic  pressure,  Pj,  is  defined  as; 


P.  =  P. 


P,gz 


+  P 


(52) 


Here,  P  '  =  total  pressure,  and 
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0,  when  i  =  l 

(Pj-p,)gh,,  when  i=2 


(53) 


Pf  = 


The  viscosity,  v„  depends  on  the  medium;  for  water,  Vj=V|,  while  for  mud,  v=v2=n*lp2  from  Eq.  21 .  The  complex 
viscosity,  n*,  depends  on  the  chosen  constitutive  model  for  mud  rheology;  see,  e.g.  Eqs.  23  through  28, 

The  following  ten  boundary  conditions  are  introduced  next: 

At  the  mud  bottom,  z=-(h|+h2) 


At  the  mud-water  interface,  z=-h,+f2 


Wj  =  0 
Uj  =  0 

Wj  =  w, 


(54) 

(55) 


(56) 

(57) 


(58) 


dw  dw 

Pi  -2p,»'|-^—  ■■P|gf2  “  P2 


(59) 


\  au, 

aw, 

f  au. 

aw,] 

1  1 

=  Pj'-a 

— :  +  — : 

az 

ax  j 

[  az 

ax  j 

At  the  free  surface,  z=J'| 


3r, 


"  at 


P.--2p...^  =  0 


(60) 


(61) 

(62) 


Under  cyclic  wave  loading,  consider: 


au, 

aw, 

dz 

ax 

=  0 


f.  =  f,(z)exp[i(kx-ait)] 


(63) 


(64) 


Uj  =  Uj{z)exp[i(kx-wt)] 
Wj  =  w,(z)expli(kx-a)t)] 
P,  =  Pj(z)exp(i(kx-u>t)] 


(65) 

(66) 
(67) 
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and  for  the  shear  stress: 


3u  flw 

_ ;  + _ ! 

dz  dx 


Pii'i 


3Ui(z) 

dz 


•Wi(z)*ik 


expli(kx-ft)t)] 


(68)  j 

I 


Implicit  solutions  for  the  wave  number  and  attenuation  coefficient,  velocity,  pressure  and  shear  stress  are 
obtained  (Dalrymple  and  Liu,  1978).  Maa  and  Mehta  (1987,1990),  Chou  (1989)  and  others  extended  this  method 
to  incorporate  stratified  mud  properties.  Jiang  (1993)  has  extended  the  approach  to  include  second  order  effects 
arising  from  finite  amplitude  waves. 


4.6  MODELING  WATER  COLUMN  SEDIMENT  DYNAMICS 

Since  we  are  focussing  on  vertical  transport  in  the  water  column,  the  basic  equation  is  the  balance  of  sediment 
mass  in  terms  of  the  settling  and  upward  mass  diffusion  fluxes,  F,  and  respectively: 

.^  +  _L(F,,-F  )=0  (69) 

at  dz^  ” 

where  Pd=4»Pj=(p-Pw)Pj/(Pj-Pw)  is  density  or  mass  concentration  of  sediment  in  suspension  (dry  sediment 

mass  per  unit  volume  of  suspension).  We  note  that  Fa=-k^dpo^dz,  where  k,  is  the  vertical  mass  diffusivity,  and 
F,=w^d*  where  w,  is  the  settling  velocity.  Given  the  comparatively  short  wave  periods  in  relation  to  the  gradually 
varied  trend  in  Po.  it  is  customary  to  consider  vertical  transport  on  a  wave-average  basis. 

The  initial  condition  for  solving  Eq.  69  is  prescribed  by  Pd(z,0),  and  at  the  water  surface,  z=0,  the  usual  no 
net  flux  condition  is,  F/0,t)-F,(0,t)=0.  Two  approaches  have  been  used  in  prescribing  the  erosion  flux  condition 
at  the  mud  surface,  z=-h,,  as  follows. 

Some  investigators,  e.g.  Thimakom  (1980;  1984)  and  Cervantes  (1987),  have  prescribed  /»D(-h,,t),  following 

a  similar  approach  commonly  used  in  coarse-grained  sediment  transport.  In  that  case  a  reference  concentration  is 

( 

prescribed  close  to  the  bed  at  elevation  z=-hi  5,,  where  5,  is  a  small  distance  above  the  bed  which  depends  on 
the  work  done  by  the  fluid  in  dislodging  and  raising  the  potential  energy  of  the  surface  particle  elevated  to  6,  (Smith, 
1977).  Since,  however,  bed  load  transport  does  not  occur  when  the  sediment  is  fine-grained,  especially  cohesive 
(Partheniades,  1977),  for  most  practical  purposes  it  is  reasonable  to  set  6,=0.  In  any  event,  defining  /3(t)=Pn(- 
h,,t)/pj,(t),  where  Pp(t)  is  the  instantaneous  depth-mean  value  of  Pod),  experimental  data  for  /3(t)  and  an  empirical- 
fit  equation  are  given  in  Fig.  20.  The  tests  (T-1 ,  T-2  and  T-3)  were  performed  in  a  wave-flume  shown  schematically 
in  Fig.  21,  using  an  estuarine  mud.  See  Section  5.2  for  a  further  description  of  the  flume.  The  wave  frequency  w, 
and  the  nominal  wave  amplitude,  Sq,  for  each  test  ire  given  in  Table  8.  Other  relevant  test  parameters  are  prescribed 

in  Cervantes  (1987).  In  the  empirical  equation  of  Fig.  20,  /8=jS/)8o,  and  t=5,t;  thus  data  fit  for  )3(t)  is 

dependent  on  0^  and  6,.  These  values  are  given  in  Table  8. 


34 


NORMALIZED  TIME,  t 


Figure  20.  Variation  of  0  with  normalized  time,  t  (after  Cervantes,  1987). 


Figure  21.  Wave-mud  interaction  flume  (after  Maa,  1986). 


Table  8.  Test  parameters  including  wave  frequency  w  and  nominal  wave  amplitude,  ao,  and  coefficients  0o.  04  and 


Test 

0) 

(rad/s) 

ao 

(m) 

00 

04 

(s') 

5. 

(s') 

T-1 

6.3 

0.030 

18.5 

0.045 

0.00055 

T-2 

3.1 

0.015 

50.0 

0.220 

0.00090 

T-3 

6.3 

0.035 

48.0 

0.060 

0.00050 

A  noteworthy  observation  relative  to  Fig.  20  is  that  0(t),  hence  0(t),  is  observed  to  rise  rapidly  initially, 
followed  by  its  f^proach  to  a  steady  state  value  close  to  one.  The  initial  rise  represents  a  mass  erosion-type  bed 
failure  phenomenon  represented  by  the  Mohr-Coulomb  equation;  see  Lamb  and  Whitman  (1969),  whose  effect  is 

mainly  embodied  in  the  coefficient  0^,  which,  for  a  given  bed,  apparently  depends  not  only  on  t^,  but  also  on  dr^dt 

(Cervantes,  1987).  In  Fig.  20,  the  trend  of  increasing  mass  erosion  with  increasing  0^  is  generally  evident.  This 
characteristic  response  of  the  bed  suggests  that  the  release  of  fine  sediment  into  the  water  column  from  the  bed  is 
likely  to  be  an  episodic  process,  with  mass-eroded  bottom  sediment  entrained  perhaps  by  coherent  bursts.  The 

approach  of  0  to  unity  suggests  that  following  the  time-dependent  effects  associated  with  experimental  start-up,  the 
vertical  dry  density  profile  became  uniform.  This  was  so  because  pD('h;.0  was  actually  measured  above  the 
lutocline.  Thus,  fluid  mud  generation  is  not  represented  in  the  character  of  0{i)  in  Fig.  20,  for  which  the  datum, 
z=-hi,  was  located  approximately  at  the  elevation  of  the  lutocline.  Thus,  in  effect  5,  was  approximately  equal  to 
the  thickness  of  the  fluid  mud  layer. 

It  can  be  easily  shown  (Cervantes,  1987)  that  by  specifying  0(t),  and  introducing  it  in  Eq.  69  renders  the 
problem  as  one  of  initial  value,  which  can  be  solved  for  pp(t)  given  pp(0),  and  it  is  recognized  that  Pp(t) 
approaches  a  constant  value,  p^j,  as  t-»oo.  The  vertical  distribution  of  Po  is  not  obtained  by  this  method. 

In  low  energy  environments  where  significant  fluid  mud  generation  usually  does  not  occur,  j3=  1  is  not  an 
unreasonable  assumption.  Elsewhere,  however,  in  general  wave-generated  fluid  mud  must  be  treated  explicitly  as 
an  eroded  sediment  mass,  especially  since  in  this  state  mud  is  relatively  easily  transported  horizontally  by  currents. 
Process  models  for  fluid  mud  generation  have  been  developed  (e.g.  Chou,  1989;  Foda  et  al.,  1993)  to  calculate  the 
fluid  mud  thickness  (see  Section  3.2). 

Erosion  of  the  fluid  mud-water  interface  occurs  as  a  process  involving  shear-induced  warping  of  the  interface, 
which  leads  to  its  destabilization  and  mixing  (Scarlatos  and  Mehta,  1993;  Adams  et  al.,  1993).  This  process  is 
mainly  dependent  on  the  Richardson  number  and  also  the  settling  velocity  of  the  particle.s  being  entrained. 
Experiments  under  uni-directional  shear  flows  have  yielded  the  following  expression; 

p,AbU  f  A  B 

e  =  -  I  —  +_ - 

g  I  Ri  Pe''^ 

where  p,  =  mean  upper  (mixed)  layer  density,  b  =  buoyancy  =  g(P2-p*)/Pw.  Pi  ~  mud  density,  p„  =  water 
density,  Ab  =  buoyancy  jump  across  the  fluid  mud-water  interface,  Ri  =  Richardson  number  =  h|Ab/U%  U  = 
mean  velocity  of  the  upper  layer,  Pe  =  Peclet  number  =  Uh|/k™,  k„,  =  molecular  diffusivity  and  A,B.D  are 
empirical  constants  (Mehta  and  Srinivas,  1993).  The  coefficient  D  depends  on  the  settling  velocity;  thus,  a  non¬ 
settling  constituent  such  as  salt  D=0,  and  Eq.  70  becomes  applicable  for  example  to  salt  entrainment.  Note  that  the 
Peclet  number  dependent  term  is  important  only  at  very  high  values  of  Ri  when  the  flow  becomes  viscous.  For  most 
practical  applications  it  is  appropriate  to  set  B=0.  The  applicability  of  Eq.  70  under  a  wave  field  remains  to  be 
demonstrated,  however. 
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The  second  approach  with  respect  to  the  erosion  flux  involves  the  specification  of  3pp(-hi,t)/3t  (Ross  and 
Mehta,  1989).  Thus,  from  Table  5, 

3pp(-h,,t)  ^ 

at  h,  Tg(t) 

in  which  Tgft)  is  obtained  from  tr(z)  via  the  continuity-based  z-t  transformation 

dz  _  app(t)  (72) 

at  Pj(z)  at 

where  pjfz)  is  the  depth-varying  bottom  mud  density.  In  the  difference  form,  Eq.  72  can  be  written  as 
At=Az/{[hi/p2(z)][app(t)/at]}.  Thus,  a  change  in  Tr  over  the  differential  scour  depth,  Az,  can  be  interpreted  in 
terms  of  the  corresponding  change  in  tr  over  At.  Note  that,  since  rg  t3q)ically  increases  with  p,  by  virtue  of  Eq.  71 
scour  is  arrested  when  the  condition,  rR=r,„  is  attained  at  a  particular  depth. 

The  datum,  z=-h,,  can  be  conveniently  located  at  the  level  between  the  bed  and  fluid  mud.  In  this  way  fluid 
mud  can  be  treated  in  an  approximate  way  as  a  suspension  formed  by  bed  erosion  specified  by  Eq.  71.  This  equation 
must  be  "calibrated"  by  adjusting  the  value  of  Mj  to  account  for  the  rate  of  generation  of  fluid  mud,  as  well  as  its 
entrainment  into  the  water  column  above  the  lutocline. 

In  high  energy  environments  wherein  a  considerable  amount  of  sediment  is  resuspended,  the  influence  of  the 
negative  buoyancy  and  turbulence  damping  characteristic  of  the  sediment-water  mixture  can  be  substantial,  and  tend 
to  retard  upward  mass  diffusion  (Wolanski  et  al.,  1988).  It  is  therefore  essential  to  amend  the  commonly  chosen 
neutral  diffusivity,  ko,  by  the  damping  function  d>  =  H-Q'5z/L,  where  aj  is  an  empirical  coefficient,  and  L  is  the 

Monin-Obukov  length  scale.  Thus,  k,=ko^,  and  a  common  empirical  form  for  is:  d'=(l+|3jRij)'“*,  where  the 
gradient  Richardson  number,  Rig=(g/p)(dp/dz)/(du/dz)^,  p  is  the  density  of  the  sediment-water  mixture  and  u(z)  is 
the  horizontal  flow  velocity.  We  therefore  have  the  diffusion  flux,  Fj=-ko(l  +/ijRi^)‘°'3pp/az.  which  is  plotted  in 

Fig.  22  against  dp^dz  for  stable  and  unstable  stratifications,  given  0^=2,  and  other  ."'evant  parameters 

based  on  the  work  of  Ross  (1988).  See,  e.g.  Horikawa  (1978)  for  a  discussion  on  the  choice  of  kj,  and  Ross  (1988) 
for  and  values. 

The  settling  velocity,  w„  in  general  varies  with  pp.  Defining  p^  such  that  for  all  Pd<Pdi.  settling  is  free,  i.e. 
the  inter-particle  collision  frequency  is  so  low  that  particles  or  aggregates  settle  essentially  independently  of  each 
other,  the  following  expressions  can  be  used: 

For  Pp<Pdi> 

w,  =  contant  (i.e.  independent  of  Pp)  (73a) 
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Figure  22.  Diffusion  flux,  F^,  as  a  function  of  dpjdz  for  and  13^=4.17  (modified  from  Ross,  1988). 


For  /Od>/}d,, 


w 


(73b) 


where  coefficients  a,  b,  m  and  n  depend  on  sediment  properties.  In  addition,  as  noted  later,  they  also  depend  on 
the  degree  of  agitation  in  the  water  column.  In  Eq.  73b,  for  relatively  low  values  of  pp,  b^>  >p^,  hence  w,  =  ab 
^>0".  This  is  called  flocculation  settling  in  which  w,  increases  with  p^  due  to  a  corresponding  increase  in  the  inter¬ 
particle  collision  frequency.  At  high  values  of  pp,  b^  <  <Po-  hence  w^  =  a/pp-" ".  In  this  case  w^  decreases  (m  >  n/2) 
with  increasing  pp  due  to  hindered  settling. 

Settling  velocity  data  that  approximately  conform  with  Eq.  73b  are  shown  in  Fig.  23  with  a=33.38,  b=2.537, 
n=1.83  and  m=1.89.  The  settling  flux,  F,,  is  obtained  from 

F  =  (4  ^8) 

’  (PD+b')™ 

The  quantity  Pm  is  somewhat  arbitrarily  prescribed  in  Fig.  23  as  0. 1  g/1,  ppj  =  2.46  g/1  is  the  value  of  pp  at  which 
the  settling  velocity  peak  occurs,  and  pp,  =  4.38  g/1  corresponds  to  the  peak  value  of  F,.  A  lutociine  occurs  when 
Pn  exceeds  Po3l=(p„-p*V(p,/Pw)l. 

For  a  given  sediment,  the  coefficients  a,  b,  m  and  n  are  found  to  be  sensitive  to  the  degree  of  agitation,  or 
turbulence,  in  water.  As  observed  from  Fig.  24,  agitation  by  vertically  oscillating  rings  in  a  laboratory  column 
caused  the  particulate  aggregates  to  rupture  and  consequently  settle  more  slowly  than  under  quiescent  conditions 
(Wolanski  et  al.,  1992). 
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Figure  23.  Settling  velocity  (and  flux)  variation  with  pp.  Data  points  obtained  from  tests  using  a  lake  (Okeechobee, 
Florida)  sediment  in  a  laboratory  column  (after  Hwang,  1989). 


Figure  24.  Settling  velocity  versus  pp  data  using  a  bay  (Townsville,  Australia)  sediment  under  two  conditions: 
quiescent  water  and  water  agitated  by  vertically  oscillating  rings  (adapted  from  Wolanski  et  al.,  1992). 
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LABORATORY/FIELD  DATA  AND  MODEL  SIMULATIONS 


5.1  INTRODUCTION 

As  noted,  given  the  compliant  and  highly  dissipative  nature  of  bottom  mud,  wave-induced  water  motion  is 
significantly  contingent  upon  the  choice  of  the  constitutive  relationship  describing  mud  rheology.  For  illustrative 
purposes  we  will  consider  mud  to  be  viscoelastic,  a  special  case  of  which  is  mud  as  a  viscous  fluid.  Note  that  mud 
motion  under  typical  conditions,  both  in  the  laboratory  and  in  the  field,  is  sufficiently  damped  by  the  high  viscosity 
of  mud,  which  tends  to  be  order  of  1(F  to  10^  times  greater  than  water.  Thus,  turbulence  within  mud  can  be  ignored 
in  most  practical  applications,  even  though  it  does  occur  to  some  degree  below  the  lutocline,  and  contributes  towards 
the  entrainment  of  fluid  mud  (Wolanski  et  al.,  1988).  Laboratory  applications  are  described  first,  followed  by  field, 

5.2  WAVE  FLUME  RESULTS 

S.2.1  Test  Facility,  Experiments  and  Modeling  Approaches 

Experiments  were  conducted  in  a  20  m  long,  46  cm  wide  and  45  cm  high  plexiglass  flume  shown  in  Fig.  21 . 
Progressive  waves  were  generated  by  a  programmable,  plunging -type  generator.  Mud  was  contained  in  a  16  cm 
deep,  8  m  long  trench  with  end  slopes  as  shown.  The  down-wave  end  of  the  flume  beach  was  artificially  roughened 
to  absorb  wave  energy  and  thus  minimize  reflection.  Details  of  flume  description  are  found  in  Ross  (1988).  Muds 
of  different  composition,  degree  of  consolidation  and  vertical  density  structure  were  tested.  Figure  25a  shows  the 
time-variation  of  stratified  density  profiles  of  a  consolidating  mud  from  Cedar  Key,  Florida  (designated  CK) 
prepared  by  depositing  the  sediment  from  an  initially  formed  suspension  in  the  flume.  Note  that,  in  general,  in  such 
experiments  the  time-scale  of  measurable  density  variation  is  considerably  greater,  on  the  order  of  tens  of  minutes, 
than  the  typical  wave  period  and,  therefore,  consolidation  does  not  measurably  contaminate  the  wave-averaging 
process  usually  inherent  in  modeling.  In  contrast  to  Fig.  25a,  Fig.  25b,  shows  a  practically  uniform  clayey  mixture 
composed  of  an  attapulgite  and  a  kaolinite  (designated  AK)  in  equal  proportions  by  weight.  The  bed  was  prepared 
by  pouring  a  pre-mixed  sediment  slurry  over  the  flume  bottom. 

All  experiments  were  conducted  with  a  nominally  constant  mud  thickness,  a  constant  water  depth,  and 
monochromatic,  progressive  waves  of  selected  frequency  and  forcing  amplitude.  Measurements  in  general  included 
density  profiles  from  the  water  surface  to  the  flume  bottom,  wave  height  decay  over  the  length  of  the  trench,  mud 
elevations  along  the  trench,  water  velocity,  mud  velocity  and  acceleration,  and  vertical  total  and  pore  pressure 
profiles.  For  the  types  of  transducers  and  measurement  techniques  used  see  Maa  (1986),  Ross  (1988)  and  Jiang 
(1993). 

The  presented  experimental  results  encompass  three  categories;  wave  attenuation,  water  and  mud  dynamics  and 
bed/interface  erosion,  in  consonance  with  Fig.  1.  With  respect  to  wave  attenuation  and  mud  dynamics  the  data  are 
compared  with  simulations  using  the  models  of  Maa  (1986)  and  Jiang  (1993).  Both  consider  water  to  be  Newtonian 
and  mud  to  be  a  viscoelastic  continuum.  The  equations  for  continuity  and  the  horizontal  and  vertical  components 
of  motion  for  these  models  are  given  by  Eqs.  49,  50  and  5 1 .  respectively.  The  basic  problem  formulation  for  both 
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Figure  25.  a)  Time-variation  of  the  bottom  mud  density  profile  during  an  erosion  experiment  (Run  1)  using 
mud  from  Cedar  Key,  Florida  (after  Maa,  1986);  b)  Uniform  density  profile  of  a  mixture  of  an  attapulgite  and  a 
kaolinite  (after  Jiang,  1993). 
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models  is  given  in  Section  4.5.2.  Differences  betw<..en  the  two  models  are  as  follows:  The  small  amplitude  solution 
of  Maa  (1986)  is  based  on  linearized  boundary  conditions.  In  Jiang  (1993),  the  non-linear,  finite  amplitude  solution 
is  carried  out  to  second  order  using  the  perturbation  expansion  technique.  In  Maa  (1986),  bottom  mud  can  be  sub¬ 
divided  into  as  many  layers  as  desired,  each  specified  by  the  density  and  the  two  viscoelastic  constants  for  the  Voigt 
constitutive  model  (Eq.  25).  Thus  stratified  mud,  such  as  shown  in  Fig.  25a,  can  be  parameterized  in  this  model. 
In  Jiang  (1993)  the  bed  properties  are  vertically  uniform,  such  as  represented  by  the  density  profile  in  Fig.  25b. 
and  the  viscoelastic  model  is  given  by  Eq.  27,  which  is  more  general,  and  approaches  the  Voigt  model  at 
frequencies  greater  than  about  10  Hz.  Neither  model  includes  interfacial  erosion;  hence  the  bottom  mud  mass  is 
conserved. 

The  viscoelastic  assumption  causes  the  kinematic  viscosity,  i<„  in  the  equations  of  motion  (Eqs.  50  and  51)  for 
mud  to  be  a  complex  number,  the  form  of  which  depends  on  the  choice  of  the  constitutive  model,  as  noted 

earlier.  In  general,  for  a  given  mud  and  temperature,  fi*  varies  with  the  forcing  amplitude,  frequency  and  mud 
density  (Jiang,  1993).  As  noted  in  Section  3.2,  in  a  ccntrolled-straii.  rheometer,  Chou  (1989)  found  the  strain 
amplitude  to  be  an  important  parameter  in  determining  mud  response  (viscous,  viscoelastic  or  elastic),  and  frequency 
to  be  of  secondary  importance  within  the  range  of  frequencies  examined.  In  a  conirolled-stress  rheometer  Jiang 
(1993)  examined  mud  response  to  oscillatory  stress  amplitude  as  well  as  forcing  frequency.  Several  muds  were 
tested,  the  stress  amplitude  was  varied  from  0.3  to  25  Pa,  and  frequency  from  0.02  to  4.5  Hz.  Typically,  for  stress 
amplitudes  ranging  from  0.6  to  10  Pa,  the  response  was  found  to  be  relatively  weakly  dependent  on  the  amplitude. 
On  the  other  hand,  frequency  dependence  was  more  significant.  As  an  example,  in  Fig.  26,  the  magnitude  of  the 
equivalent  viscosity,  |/i*|  =  IpjI'*!  =P2(p'-+p''-)'”,  is  plotted  againsi  frequency,  f,  for  a  mud  from  offshore  of 
Mobile  Bay,  Alabama  (designated  MB),  the  site  of  a  dredged  material  deposit  from  the  bay  (Mehta  and  Jiang, 
1993).  The  strong  dependence  of  |  p*  |  on  f  is  believed  to  be  due  to  the  thermodynamic  process  associated  with  the 
vibratory  response  of  mud  to  cyclic  forcing;  at  low  frequencies  the  process  is  probably  isothermal,  while  at  high 
frequencies  it  tends  to  become  adiabatic  (Krizek,  1971;  Schreuder  ei  al.,  1986).  This  change  evidently  has  the 
overall  effect  of  decreasing  |p*j  with  increasing  f.  Increasing  <p  in  turn  increases  the  stiffne.ss  and  the  effective 
viscosity. 

The  effective  viscosity  magnitude  in  Fig.  26  has  been  determined  for  the  viscoelastic  model  shown  in  Fig.  18a, 
and  is  thus  defined  by  p,  G,  and  Gj.  These  three  parameters  were  experimentally  found  to  depend  on  f  according 
to: 

p,  G,  or  Gj  =  exp(f)-f'^  05) 

where  values  of  the  coefficients  e  and  A  are  given  in  Table  9  (Mehta  and  Jiang.  1993).  Relatioi.'hips  such  as  the 
ones  in  Fig.  26  inherently  suggest  that,  for  example,  as  waves  pass  over  mudbanks  the  wave  spectral  .hape  may 
be  drastically  altered  due  to  frequency-selective  attenuation,  as  observed  for  example  off  Surinam  in  Fig.  27  (Wells 
and  Kemp,  1986). 
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Figure  26.  Magnitude  of  equivalent  dynamic  viscosity  (applicable  to  the  viscoelastic  model  of  Fig.  18a)  as  a 
function  of  wave  forcing  frequency  and  mud  volume  fraction,  <(>,  for  MB  mud. 

Table  9.  Coefficients  of  Eq.  5.1  for  MB  mud 
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-0.030 
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-0.696 

8.374  , 
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Wave-induced  erosion  is  separately  examined  in  Section  5.2.3  via  the  approaches  described  in  Section  4.6, 
following  wave  attenuation  and  associated  mud  dynamics.  The  first  approach  yields  the  time-variation  of  the  depth- 
averaged  suspended  sediment  concentration  or  dry  density,  during  erosion,  while  the  second  gives  the  time- 
variation  of  the  vertical  profile  of  concentration,  Pi,(z,t).  Some  relevant  laboratory  results  are  considered  next. 

5.2.2  Wave  Attenuaiion  and  Mud  Motion 

In  Fig.  28  the  wave  attenuation  coefficient,  k,=-(l/x)ln(a,/ao),  measured  over  a  distance  approximately  equal 
to  the  length  of  the  mud  trench  (Fig.  21),  has  been  plotted  against  wave  frequency,  to,  for  AK  mud.  The  depth  of 


Figure  27.  Wave  spectra  from  synchronous  wave  stations  22,  1 1  and  4  km  offshore,  showing  Joss  of  wave  energy 
and  change  in  the  shape  of  the  spectrum  over  a  fluid-like  mud  bottom  along  the  central  Surinam  coast.  =  root 
mean  square  wave  height,  T  =  wave  period  and  h  =  water  depth  (after  Wells  and  Kemp,  1986). 

water  (h,)  in  these  tests  was  16  cm,  while  the  mud  depth  (h,)  was  varied.  The  input  wave  amplitude  {=a„)  was 
1.3  cm.  In  addition  to  the  effect  of  water  depth  on  the  degree  of  attenuation,  the  maximum  in  the  k,-w  curve  is  an 
interesting  and  characteristic  feature  that  requires  an  explanation.  In  the  special  case  of  shallow  water  and  viscous 
mud  (see  Section  4.5),  Jiang  and  Mehta  (1992)  showed  that  the  degree  of  attenuation,  as  reflected  by  k,,  depends 
on  the  ratio,  hj/S,,,,  where  6„=(2v2/w)''^  represents  the  viscous  mud  boundary  layer  thickness.  The  depth-mean  value 
of  the  velocity  gradient,  dulBz,  and  in  turn  the  rate  of  wave  energy  dissipation  in  mud,  is  maximum  when  h;/6„,  is 
equal  to  or  is  close  to  one,  and  decreases  as  hj/S,,,  increases  above  one  or  decreases  below  one.  This  criterion  is 
modulated  when  mud  rigidity  is  measurable.  Thus,  MaePherson  (1980)  showed  that  the  peak  wave 
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Figure  28.  Wave  attenuation  coefficient  against  frequency  for  two  AK  mud  =  0.16)  depths,  12  cm  and  18 
cm.  Water  depth  was  16  cm.  The  input  wave  amplitude  was  1.3  cm.  Circles  are  data;  lines  are  simulations 
based  on  the  model  of  Jiang  (1993). 


ACCELERATION  AMPLITUDE  (m/s2) 

Figure  29.  Comparison  between  measured  and  model-calculated  profiles  of  the  amplitude  of  horizontal  acceleration 
in  a  water-mud  (AK)  system.  The  water  depth  was  16  cm,  mud  thickness  18  cm,  water  surface  forcing  wave 
amplitude  0.5  cm,  wave  frequency  1  Hz  and  </>  =  0.12  (after  Jiang,  1993). 
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attenuation  rate  occurs  when  v  =  «'2/(ghi’)*'^  and  G  =G/p2ghi  are  of  the  same  order  of  magnitude.  In  any  event  it 
is  evident  and  noteworthy  that,  as  a  result  of  the  strong  dependence  of  k*  on  w,  the  mode  of  change  in  the  wave 
spectral  shape  as  the  waves  pass  over  a  mudbank  will  critically  depend  on  the  character  of  the  input  spectrum,  mud 
viscosity  and  mud  depth. 

Based  on  the  results  from  another  test  using  the  same  mud  (AK),  a  comparison  between  the  measured  and 
model-calculated  amplitude  of  the  horizontal  acceleration.  Uj,  within  mud  is  shown  in  Fig.  29.  Note  the  drastic 
reduction  in  mud  acceleration  relative  to  that  in  the  water  column  due  to  the  dissipative  nature  of  the  bottom 
medium. 

A  noteworthy  feature  of  bottom  mud  subjected  to  non-linear  waves  is  residual  mud  mass  transport,  which  in 
a  large  number  of  situations  is  a  significant  factor  leading  to  sedimentation,  e.g.  in  port  channels  and  basins.  The 
Lagrangian  mean  particle  drift,  or  Stokes’  drift,  Ul,  associated  with  "non-closed"  trajectories  of  mud  particles  must 
therefore  be  calculated.  In  general,  the  velocity  of  a  particular  particle  with  a  mean  position  of  (x,,z,)  is  u(x,-t-f ', 
3^1 +£').  where  f '  and  are  instantaneous  position  coordinates  along  the  particle  trajectory.  Based  on  Taylor  series 
expansion,  we  have  (Sakakiyama  and  Bijker,  1989;  Jiang,  1993); 

du,  f  3u, 

UL(X,+f',Z,^r)  =  ^  *  -gj 

where  overbars  denote  time-averaging  over  the  wave  period  anduj  is  the  Eulerian  mass  transport  velocity.  To  obtain 
Ul  theoretically,  the  wave-averaged  intergals  and  derivatives  on  the  right  hand  side  of  Eq.  76  must  be  obtained 
through  an  appropriate  hydrodynamic  model.  In  Fig.  30a  comparison  is  made  between  measurement  and  prediction 
for  AK  mud,  using  the  model  of  Jiang  (1993).  The  measurements  were  made  by  observing  the  horizontal 
displacement  of  an  initially  vertically  injected  dye  streak  adjacent  to  the  flume  side  wall  made  of  plexiglass.  Two 
sources  of  error  are  noteworthy  in  terms  of  the  observed  discrepancy  between  the  data  points  and  the  model- 
calculated  profile.  Firstly,  the  data  points  incorporate  the  effects  of  drag  close  to  the  side  wall  and,  secondly,  the 
calculated  mass  transport  velocity  profile  is  approximate  due  to  analytical  limitations  (Jiang,  ^993).  Nevertheless, 
the  overall  trends  do  compare  favorably.  Accurate  field  information  on  mass  transport  velocity  is  apparently  sparse; 
however,  the  model  was  applied  to  predict  the  rate  of  mass  transport  of  transient  mudbanks  (of  KI  mud)  off  the 
coast  of  Kerala  (India).  The  predicted  0.22  km/d  shoreward  transport  rate  at  the  onset  of  monsoon  at  least  appears 
to  be  of  the  correct  order  of  magnitude  (Mathew,  1992;  Jiang,  1993). 

From  an  experiment  meant  to  investigate  mud  dynamics,  measured  and  predicted  profiles  of  velocity  and 
dynamic  pressure  are  compared  in  Fig.  31a,b  for  CK  mud.  The  calculations  are  based  on  the  model  of  Maa  (1986), 
in  which  the  mud  bottom  was  subdivided  into  four  layers  of  thicknesses  indicated,  each  characterized  by  density, 
viscosity  and  shear  modulus  of  elasticity.  The  dynamic  pressure  at  first  decreases,  then  increases  with  depth  below 
the  interface.  This  variation  has  been  explained  by  Maa  (1986)  in  terms  of  the  relative  influences  of  the  effective 
viscous  term  and  the  pressure  term  in  the  equation  of  motion.  Note  that  the  model  predicts  an  interfacial  shear  stress 
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Figure  30.  Mass  transport  velocity  in  AK  mud.  Wave  amplitude  2  cm,  frequency  1  Hz,  h,  =  14  cm  and  mud 
thickness  =  17  cm  (after  Jiang,  1993). 


AMPLITUDE,  (p/piga) 


Figure  31.  a)  Horizontal  and  vertical  velocity  amplitudes  in  water  and  CK  mud;  b)  Dynamic  pressure  amplitude 
profile  corresponding  to  Fig.  31a.  Smooth  profiles  have  been  drawn  based  on  spatially  discretized  calculations 
(after  Maa  and  Mehta,  1987;  Maa,  1986). 
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(amplitude,  t^)  of  0.19  Pa,  which  was  too  small  to  cause  a  substantial  amount  of  interfacial  erosion  during  the 
experiment.  Yet  the  bottom  mud  showed  a  measurable  horizontal  oscillation  driven  by  the  pressure  gradient.  This 
type  of  a  situation  is  likely  to  be  ubiquitous  in  nature,  with  the  result  that  under  typically  gentle,  fair  weather  wave 
action  mud  may  oscillate,  even  while  the  upper  water  column  is  relatively  free  of  turbidity,  since  the  interfacial 
stress  is  insufficient  to  cause  significant  mud  scour.  Thus,  for  example,  in  Lake  Okeechobee,  Florida,  daily  wave- 
induced  mud  oscillation  tends  to  retain  the  top  5-20  cm  layer  of  bottom  mud  in  a  fluid-like  state,  which  in  turn 
presumably  facilitates  the  transport  of  nutrients  across  the  mud-water  interface,  and  influences  the  trophic  state  of 
this  large  and  shallow  lake  (Jiang  and  Mehta,  1992). 

5.2.3  Bed/Interface  Erosion 

The  time-variation  of  the  wave-mean  and  depth-mean  dry  density  (concentration)  of  suspended  sediment  of  a 
bed  of  CK  mud  is  shown  in  Fig.  32  in  test  T-1  of  Cervantes  (1987)  mentioned  in  Section  4.6.  The  normalized 
parameters  are:  dry  density,  Po^’^ds-  ^=(w5/h|)t.  Note  the  almost  instantaneous  rise  of  the 

normalized  dry  density  to  0.86,  followed  by  a  drop  and  subsequent  attainment  of  a  steady  state  (p„-*l ) .  The  peculiar 

nature  of  the  curve  reflects  the  corresponding  variation  of  /3  with  normalized  time  (t  )  shown  in  Fig.  20.  The 
function  given  in  the  inset  of  Fig.  20  has  two  components,  the  first  of  which  reflects  surface  erosion  and  the  second 
is  attributed  to  Mohr-Coulomb  type  bed  failure  or  mass  erosion.  Initially  the  influence  of  the  second  term  is 
dominant  in  Fig.  32  and,  as  steady  state  approaches  surface  erosion  (first  term)  becomes  the  main  contributor.  Thus, 
qualitatively  the  response  curve,  i.e,  the  time-variation  of  p^,  can  be  divided  into  the  three  phases  shown  in  Fig.  33. 
In  phase  P-1  mass  erosion  is  dominant,  P-II  is  characterized  by  settling,  and  surface  erosion  is  dominant  in  P-lII. 
The  occurrence  of  the  mass  erosion  peak  may  be  of  noteworthy  significance  in  nature.  Although  supportive  field 
data  are  apparently  sparse,  such  peaks  are  probably  quite  common  at  the  time  of  comparatively  sudden  onset  of 
increased  storm  wave  action,  amounting  to  an  episodic  generation  of  turbidity. 

The  time-variation  of  the  vertical  profile  of  the  dry  density  during  an  erosion  experiment  is  shown  in  Fig.  34. 
The  bed  was  composed  of  HB  (from  Hillsboro  Bay,  Florida)  mud,  and  simulation  is  via  the  m6del  of  Ross  (1988), 
which  is  based  on  the  sediment  transport  equation  given  in  Section  4.6.  Note  the  rapid  generation  of  the  lutocline, 
similar  to  what  occurs  in  the  field,  e.g.  Fig.  35  (adapted  from  Kemp  and  Wells,  1987).  Of  the  four  instantaneous, 
vertical  dry  density  profiles  for  suspended  sediment  shown  in  Fig.  35,  one  represents  pre-frontal  wave  conditions, 
two  were  obtained  during  frontal  passage,  and  the  last  is  post-frontal.  The  data  were  obtained  at  a  site  on  the  eastern 
margin  of  the  Louisiana  chemier  plain  where  the  tidal  range  is  less  than  0.5  m  (Kemp  and  Wells,  1987).  Wave 
height  during  the  winter  cold  front  passage  was  on  the  order  of  13  cm  and  period  7  s.  Of  particular  interest  is  the 
development  of  a  fluid  mud  layer  by  the  frontal  wind-generated  waves  (profiles  2  and  3),  which  was  previously 
absent  (profile  1).  The  post-frontal  profile  4  further  suggests  that  this  layer  may  have  persisted  following  front 
passage,  conceivably  due  to  the  comparatively  low  rate  at  which  the  layer  was  dewatered.  The  suspension  density 
in  the  upper  water  column  was  higher  following  the  front  than  that  during  the  front,  possibly  due  to  advection  of 
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Figure  32.  Time-variation  of  normalized  dry  density  (concentration),  PdCO  =  Pt,(t)/pBs ,  forCK  mud,  during  test  T-1. 
Mud  depth  (hj)  was  16  cm,  water  depth  (h,)  was  17  cm,  forcing  wave  amplitude  (ao)  was  3  cm  and  frequency  (0 
was  1  Hz  (after  Cervantes,  1987). 


Figure  33.  Phases  in  the  time-variation  of  p^ft),  such  as  in  Fig.  32  (after  Cervantes,  1987). 
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Figure  34.  Measured  and  model-calculated  dry  density  (concentration)  profiles  during  erosion  of  Hillsboro  Bay  mud 
(HB).  Water  depth  (hi)  31  cm,  mud  depth  (hj)  12  cm,  forcing  wave  amplitude,  a<,  =  3  cm  and  wave  frequency,  f 
=  1  Hz  (after  Ross,  1988). 


Figure  35.  Vertical  suspended  sediment  dry  density  (concentration)  profiles  obtained  before,  during,  and  after 
the  passage  of  a  winter  cold  front  at  a  coastal  site  in  Louisiana.  Times  are  relative  to  time  of  measurement  of  the 
pre-frontal  profile  1  (adapted  from  Kemp  and  Wells,  1987). 
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suspended  sediment  from  a  neighboring  area  of  high  turbidity.  Also  note  that  in  Fig.  34  as  well  as  Fig.  35,  the  lack 
of  significant  upward  transport  of  sediment  above  the  lutocline  is  consistent  with  the  typically  low  wave-induced 
mass  diffiisivities  (Sheng,  1986). 

5.3  LAKE  OKEECHOBEE,  FLORIDA 
5.3.1  Setting 

A  large  part  of  the  central  basin  of  Lake  Okeechobee  in  southcentral  Florida  is  overlain  by  a  thin  layer  of  dark, 
organic-rich  (about  40%  by  weight)  mud  which  is  believed  to  be  a  significant  storage  for  externally  loaded 
phosphorus.  Figure  36  shows  the  mud  thickness  to  be  80  cm  at  most,  generally  occupying  the  deeper  part  of  the 
lake.  A  typical  variation  of  mud  density  with  depth  is  shown  in  Fig.  37.  As  observed  the  top  -  10  cm  of  the  mud 
layer  was  comparatively  low  in  density.  Below  this  layer  a  rapid  increase  in  density  occurred  and  led  to  a  relatively 
uniformly  dense  mud  underneath.  The  low  density  structure  at  the  top  is  apparently  characteristic  of  the  organic-rich 
mud  of  comparatively  uniform,  fluid-like  consistency  (Hwang,  1989). 


Figure  36.  Mud  thickness  contour  map  of  Lake  Okeechobee,  Florida  (after  Kirby  et  al.,  1989). 
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Figure  37.  Typical  mud  bottom  density  profile  based  on  vibrocore  data  from  two  nearby  sampling  sites  (after 
Hwang,  1989). 

5.3.2  Mud  Motion 

In  order  to  examine  the  nature  of  mud  motion  induced  by  wave  action  in  the  lake,  an  instrumented  tower  was 
deployed  at  a  site  near  Port  Mayaca  (Fig.  36).  The  mean  water  depth  was  1 .43  m  and  mud  thickness  0.55  m.  Water 
surface  variation  was  measured  with  a  pressure  transducer,  and  an  electromagnetic  current  meter  was  used  to 
measure  wave-induced  water  motion  87  cm  below  the  mean  water  surface.  To  measure  mud  motion  a  small 
accelerometer  was  embedded  20  cm  below  the  mud  surface  where  the  density  was  1,180  kg/m’.  The  shallow  water 
analytical  model  of  Jiang  (1993)  described  in  Section  4.5.1  was  used  to  aid  in  interpreting  the^data.  In  that  model, 
the  water  column  is  considered  to  be  inviscid,  and  mud  is  assumed  to  be  a  highly  viscous  fluid.  Mud  viscosity  was 
measured  in  a  rheometer;  the  kinematic  viscosity  at  the  density  of  1,180  kg/m’  was  found  to  be  1.76  x  10  ’  m’/s, 
which  is  very  large  compared  to  water  (10  *  mVs). 

Figure  38a  shows  a  typical  wave  spectrum  obtained  under  moderate  breeze  ( -  20  km/hr)  from  the  westerly 
direction.  The  dominant  wave  frequency  was  0.42  Hz,  with  a  significant  height  of  8  cm.  It  is  noteworthy  that 
calculations  (Mehta  and  Jiang,  1990)  showed  that  at  20  km/hr  wind  speed,  the  wave  height  would  have  been  four¬ 
fold  larger  over  an  otherwise  comparable  rigid  bottom. 

In  Fig.  38b  the  measured  horizontal  water  velocity  amplitude  spectrum  is  compared  with  that  simulated  by  the 
model  using  the  wave  spectrum  of  Fig.  38a  as  input.  Given  the  model  limitations,  particularly  with  respect  to  the 
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Figure  38.  Data  and  model  application.  Lake  Okeechobee,  Florida;  a)  Measured  water  wave  spectrum;  b) 
Comparison  between  measured  and  simulated  water  velocity  spectra;  c)  Comparison  between  measured  and 
simulated  mud  acceleration  spectra  (after  Jiang,  1993). 

assumption  of  mud  as  a  purely  viscous  fluid,  and  errors  inherent  in  the  data  collection  effort  (Mehta  and  Jiang. 
1990),  the  agreement  between  data  and  prediction  appears  to  be  acceptable. 

Referring  to  Fig.  38c  in  which  measured  and  simulated  horizontal  acceleration  spectra  in  mud  are  compared, 
the  data  as  well  as  model  simulation  indicate  the  occurrence  of  measurable  acceleration  at  about  0.04  Hz.  This 
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behavior  can  be  attributed  to  a  forced  second  order  wave  that  results  from  an  interaction  between  different 
frequencies  of  the  forcing  wind  wave,  as  for  example  occurs  in  the  case  of  surf  beat  over  relatively  flat  open  coast 
beaches  (Jiang  and  Mehta,  1992).  Note  that  dominant  seiching  in  Lake  Okeechobee  occurs  at  a  much  lower 
frequency,  on  the  order  of  10^  Hz.  While  energy  contribution  at  0.04  Hz  frequency  is  also  found  in  the  wave 
spectrum  (Fig.  38a)  and  in  the  water  velocity  spectrum  (Fig.  38b),  the  mud  acceleration  spectrum  is  enhanced  at 
the  0.04  Hz.  in  comparison  with  the  spectrum  at  the  0.42  Hz  forcing  frequency.  This  relative  enhancement  is  due 
to  the  strong  frequency  dependence  of  the  wave  attenuation  coefficient,  k,.  Model  results  for  example  yielded  k,  = 
0.0034  m '  at  0.42  Hz  and  0.0013  m  ‘  at  0.04  Hz.  Note  that  if  the  bottom  were  rigid,  k,  would  be  0(10^)m ',  which 
is  considerably  smaller. 

5.3.3  Resuspension 

Using  the  vertical  sediment  mass  transport  numerical  model  of  Ross  (1988),  the  time-evolution  of  the  vertical 
dry  density  (concentration)  piofile  was  simulated  under  selected  storm  wave  conditions.  Profiles  such  as  the  one 
in  Fig.  37  were  used  to  represent  bed  density  variation  with  depth.  For  a  given  wave  period  and  height,  the  bottom 
stress  amplitude  was  calculated  from:  r|,=0.5  f„/5|U(-h„t)^,  where  f„  is  the  wave  friction  factor  (Jonsson,  1966).  The 
bottom  velocity,  u(-h,,t),  was  obtained  from  the  linear  wave  theory  using  wave  amplitude,  frequency  and  water 
depth  as  input  parameters.  The  neutral  vertical  mass  diffusivity,  Iq,,  (Section  4.6)  was  obtained  from 
kfl=a'a^w(sinh^k(z+h|)/2sinh^kh,},  based  on  the  work  of  Hwang  and  Wang  (1982),  where  a'  is  a  free  coefficient. 
The  ko  value  was  weighted  by  the  damping  parameter  $  (Section  4.6).  Settling  velocity  data  are  given  in  Fig,  23. 

In  Fig.  39a  the  manner  in  which  the  suspension  profile  would  evolve  over  a  48  hr  period  starting  with  a  clear 
water  column  is  exemplified  for  a  6  s  (w  =  1.1  rad/s)  storm  wave  of  1  m  height  in  a  4.6  m  deep  water  column 
characteristic  of  the  central  part  of  the  lake.  Note  the  rather  rapid  formation  of  a  fluid  mud  layer  represented  by 
the  lutocline.  The  layer  thickness  was  found  to  be  about  7  cm  at  48  hr  (Hwang,  1989).  Note  also  however,  that  the 
corresponding  bed  scour  depth  can  be  shown  to  be  on  the  order  of  a  centimeter  only.  Thus  a  very  thin  bed  layer 
generated  an  order  of  magnitude  thicker  fluid  mud  layer.  Such  layers,  characterized  by  a  significant  density  gradient, 

i 

have  been  observed  in  wave  dominated  environments  (e.g.  Fig.  35).  Observe  in  Fig.  39a  the  very  slow  rise  in  the 
suspension  density  in  the  upper  water  column  due  to  characteristically  low  rates  of  upward  mass  diffusion.  At  48 
hr  the  dry  density  at  the  surface  rose  to  only  about  100  rag/1,  very  close  to  102  mg/1  reported  by  Gleason  and  Stone 
(1975)  in  the  southern  part  of  the  lake  during  a  storm.  If  wave  action  were  discontinued  at  48  hr,  the  profiles  of 
Fig.  39b  would  result  due  to  sediment  settling  under  calm  conditions.  At  72  hr,  or  24  hr  after  waves  ceased,  most 
of  the  material  is  observed  to  have  settled  to  the  bottom,  and  the  water  column  clarified. 

5.4  SOUTHWEST  COAST  OF  INDIA 

An  example  of  a  drastic  reduction  in  the  wave  energy  due  to  wave  passage  over  seasonal  coastal  mudbanks  is 
briefly  noted  here.  Figure  40  shows  the  region  offshore  of  the  town  of  .MIeppey  in  State  of  Kerala  on  the  southwest 
coast  of  India,  where  wave  spectra  were  obtained  at  the  two  sites  shown  (Mathew,  1992).  A  schematic  profile  of 
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Figure  39.  a)  Simulated  time-evolution  of  suspension  profile  due  to  0.9  m  high,  4  s  waves  in  a  4.6  m  deep 
water  column;  b)  Settling  of  sediment  once  waves  cease  1 1  hr  from  start  of  wave  action  (after  Hwang,  1989). 

the  mudbank  is  shown  in  Fig.  41  (Nair,  1988),  and  the  rheological  characteri.stics  of  the  mud  (Kl)  are  given  in  Fig. 
26  (and  Table  9).  The  bottom  beneath  the  mud  elsewhere  in  this  area  is  sandy.  Thus  the  subaqueous  bottom  mud 
is  not  always  visible  from  the  sandy  beachfacc.  even  though  the  leading  edge  of  the  mudbank  tends  to  occur  quite 
close  to  the  shoreline.  The  mean  depth  of  water  at  the  offshore  site  was  10  m.  and  at  the  inshore  site  a  little  over 
5  m.  In  fair  weather  the  mudbank  was  well  seaward  of  the  offshore  site,  in  deeper  water.  During  the  monsoon  the 
nearshore  edge  of  the  bank  was  shoreward  of  the  inshore  site,  while  the  seaward  edge  of  the  bank  was  close  to  the 
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Figure  40.  Coastal  site  off  Alleppey  in  Kerala,  India,  where  monsoonal,  mudbanks  occur.  The  pier  is  300  m 
long  (after  Mathew,  1992). 

Pertpherai  Wave 


Figure  41.  Schematic  profile  of  mudbank  off  the  coast  of  Kerala,  India  (after  Nair,  1988). 

offshore  site.  Thus  the  bottom  at  the  offshore  site  was  devoid  of  mud.  Figu.v-s  42a, b  respectively  show  examples 
of  wave  spectra  in  the  absence  and  in  the  presence  of  1  m  thick  mudbank  at  the  inshore  site.  In  fair  weather  the 
energy  reduction  is  observed  to  have  been  negligible.  (Note  the  slight  phase  shift  of  the  inshore  spectrum  relative 
to  offshore.  This  shift  may  be  an  artifact  of  the  data  analysis  procedure.)  In  the  monsoon  case  the  energy  reduction 
between  the  two  stations  was  about  85%.  In  fact,  as  a  result  of  the  mudbank  the  wave  energy  at  the  inshore  site 
was  lower  in  monsoon  than  in  fair  weather,  despite  the  occurrence  of  considerably  more  inclement  offshore  wave 
activity  during  the  monsoon.  Not  surprisingly,  monsoonal  mudbanks  in  Kerala  have  served  as  open  coast  havens 
for  small  fishing  vessels  that  would  otherwise  require  sheltered  harbors  (Nair.  1988) 
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Figure  42.  Offshore  and  inshore  wave  spectra  off  Alleppey;  a)  without  mudbank;  b)  with  mudbank  (after 
Mathew,  1992). 


The  model  of  Jiang  (1993)  was  applied  to  the  Alleppey  data  on  a  frequency -by -frequency  basis  by  selecting 
0.0025  mVs  for  the  eddy  diffusivity, for  the  water  column,  a  typical  value  (Lick,  1982),  and  taking  the  offshore 
wave  spectrum  as  input  at  x=0.  Note  that  in  Eqs.  50  and  5J,  /',  =  y,  is  replaced  by  j/i-t-f,  to  account  for  vertical 
diffusion  due  to  turbulence  in  the  water  column  in  a  very  approximate  way.  The  fair  weather  case  (Fig.  42a)  was 
treated  by  selecting  arbitrarily  high  values  of  the  rigidities  G,  and  G,  (applicable  to  Eq.  27  chosen  as  the  constitutive 
model  for  mud  rheology),  in  order  to  simulate  a  rigid  (sandy)  bottom.  For  the  (1  m  thick)  mud  bottom  case,  values 
of  fijiOf  G|(f)  and  G2(f)  for  KI  mud  (0=0. 12)  were  inputted  (Mehta  and  Jiang,  1993).  Results  shown  in  Figs.  43a, b 
indicate  that  the  inshore  wave  spectra  could  be  simulated  reasonably  well  for  the  fair  weather  (Fig.  43a)  as  well  as 
monsoon  (Fig.  43b)  scenarios.  As  noted  previously  the  model  yielded  a  depth-mean  residual  mud  velocity  (Stokes’ 
drift)  of  0.22  km/day,  which  seems  to  be  consistent  with  the  time-scale  of  mud  motion  at  the  onset  of  monsoon, 
even  though  no  detailed  mud  velocity  measurement  are  available  for  comparison  purposes  (Mathew,  1992). 

5.5  MOBILE  BERM,  ALABAMA 

The  sediment  placement  site  for  the  Mobile  underwater  berm  off  Dauphin  Island  in  Alabama,  created  as  a 
national  demonstration  project  to  highlight  the  beneficial  role  (wave  height  reduction  in  this  case)  of  dredged 
material,  is  shown  in  Fig.  44.  The  material  was  derived  from  the  ship  channel  within  Mobile  Bay  (McLellan  et  al., 
1990;  Hands,  1990).  Examples  of  offshore/inshore  wave  spectra  are  given  in  Figs.  45a,b  under  two  different 
offshore  wave  conditions.  Wave  energy  reductions  were  significant,  29%  and  46%  respectively.  Using  a  fixed-bed 
hydrodynamic  model  McLellan  et  al.  (1990)  showed  that  assuming  a  rigid  (i.e.  fixed  bed)  berm  crest  produced 
negligible  wave  energy  dissipation.  It  is  believed  that  as  a  result  of  the  shallow  depth  of  water,  on  the  order  of 


57 


Figure  43.  Comparison  between  measured  and  model-simulated  inshore  wave  spectra  off  Alleppey  in  Kerala. 
India:  a)  fair  weather  condition  (mudbank  absent);  b)  monsoonal  condition  (mudbank  present). 

6.5  m  over  the  berm,  the  main  cause  of  damping  is  energy  absorption  by  the  deposited  mud.^  Diver  observations 
at  the  site  suggested  the  occurrence  of  surface  wave-forced  interfacial  mud  waves  propagating  along  the  compliant 
crest  (McLeilan,  personal  communication).  Such  a  movement  can  be  construed  as  a  manifestation  of  the  participation 
of  the  bottom  material  in  the  energy  dissipation  process.  Furthermore,  the  high  degree  of  stability  of  the  berm 
against  displacement  and  deterioration  pointed  to  the  fact  that  wave  action  became  sufficiently  weak  over  the  berm 
to  prevent  significant  scour  of  the  berm  (Dredging  Research  Technical  Notes,  1992). 

The  wave-mud  interaction  model  of  Jiang  (1993)  was  used  to  simulate  the  inshore  wave  spectra  at  the  Mobile 
berm  site,  given  measured  offshore  spectra  in  Figs.  45a,b.  The  computations  were  done  on  a  frequency-by- 
frcquency  basis.  A  6  m  berm  height  ( =hj)  and  the  6.5  m  water  depth  ( =  h,)  were  selected  as  representative  vertical 
dimensions.  In  general,  the  effective  berm  width  for  calculating  the  degree  of  wave  attenuation  would  depend  on 
the  direction  of  wave  approach;  for  the  present  purpose  a  width  of  140  m  was  selected,  giving  some  allowance  for 
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Figure  44.  Construction  site  {2,750  m  long  corridor  for  dredged  material  placement)  for  the  Mobile  berm,  and 
offshore/inshore  sites  of  wave  measurement  (after  McLellan  et  al.,  1990). 

the  berm  slopes.  Thus  dissipation  was  assumed  to  occur  over  this  distance  only.  Rheological  parameters  for  the  mud 
(MB)  are  given  in  Table  9  (for  <i)=0.11,  corresponding  to  the  in  situ  mud  density). 

Comparisons  between  simulated  and  measured  inshore  spectra  are  given  in  Figs.  46a, b.  The  simulated  spectral 
energy  is  observed  to  be  generally  in  agreement  with  the  measurement  at  low  frequencies,  but  is  consistently  lower 
than  measured  for  frequencies  exceeding  about  0.25  Hz,  i.e.  periods  smaller  than  4  s.  In  that  context  it  must  be 
noted  that,  as  a  rule,  the  output  spectrum  from  the  model  will  have  the  same  general  shape  as  the  offshore  one, 
which  also  is  the  case  here.  Compare,  for  example,  the  measured  offshore  spectrum  with  the  simulated  inshore  one 
in  Fig.  46a.  These  two  are  generally  similar  in  shape.  Next  compare  the  measured  inshore  spectrum  with  the 
measured  offshore  one  in  the  same  figure.  Note  that  for  frequencies  greater  than  about  0.30  Hz  the  measured 
inshore  wave  energy  was  actually  slightly  higher  than  that  at  the  offshore  location.  This  feature  suggests  that  the 
source  of  the  high  frequency  inshore  waves  may  have  been  at  least  partly  different  from  that  represented  by  the 
measured  offshore  spectrum.  In  any  event,  diagnostic  applications  such  as  this  one  provide  a  basis  for  developing 
design  guidelines  for  wave  absorbing  berms  in  future.  A  hypothetical  design  example  has  been  discussed  by  Mehta 
and  Jiang  (1993). 
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Figure  45.  Offshore  and  inshore  wave  spectra  at  the  Mobile  berm  site  for  two  different  wave  conditions  at  the 
offshore  site  characterized  by  the  maximum  wave  height.  (a)  =  0.9  m,  and  (b)  =  1.5  m  (after 

McLellan  et  al.,  1990). 
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Figure  46.  Comparison  between  measured  and  model -simulated  inshore  wave  spectra  at  the  Mobile  berm  site 
corresponding  to  two  offshore  wave  conditions;  (a)  =  0.9  m  (see  Fig.  45a),  and  (b)  =  1.5  m  (see  Fig. 

45b).  Measured  offshore  spectra  have  been  included  for  reference. 
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APPENDIX  A;  EXPLANATORY  NOTES  FOR  TABLE  7 


Note  I:  The  solution  is  based  on  a  two-fluid  system  consisting  of  an  inviscid  fluid  overlying  a  lower  viscous  mud 
layer  in  laminar  flow,  which  is  assumed  to  be  infinite  in  depth,  driven  by  linear  water  waves,  A  boundary  layer 
approach  via  prescription  of  a  boundary  layer  correction  to  the  potential  in  the  mud  region  is  then  adopted 

whereby  mud  flow  is  assumed  inviscid  except  in  the  boundary  layer  region.  Two  wave  modes  emerge  from  the 
solution.  The  first  mode,  u)^=gk,,  which  is  only  true  for  infinitely  deep  mud  layer,  precludes  a  discontinuity  of 
horizontal  velocity  across  the  interface  and,  hence,  there  is  no  boundary  layer  formation  and  no  associated  damptng. 
For  the  second  mode,  ci)Vgk,  =  [0o2//3,)-lli2mhk,h,/((/)2/p,)+tanhlc,h||,  the  wave-mean  rate  of  viscous  energy 
dissipation  per  unit  area  in  the  boundary  layer  is  computed  from: 
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The  relative  dominance  of  the  two  wave  modes  is  decided  by  the  nature  of  wave  generation  (Dean  and  Dalrymple. 
1984).  While  the  first  mode  i«  more  likely  under  the  circumstance  of  waves  propagating  into  a  muddy  region,  the 
second  mode  will  dominate  if  the  waves  are  generated  at  the  interface  by  a  displacement  of  the  mud. 


Note  2:  The  solution  is  based  on  .  i  vo-lavered  viscous  fluid  system  driven  by  linear  water  waves  propagatUig  in 
water  of  constant  depth,  h,,  overlying  a  viscous  mud  layer  of  thickness,  with  appropriate  boundary'  layer 
approximation  to  render  the  solution  analytical .  These  .solutions  apply  when  the  parameters,  r,  =  («'/.(' h vjui)'" ,  i  =  1 . 
2,  are  small  whence  the  motion  is  essentially  irrotational  except  near  boundaries  where  the  vi  :ous  boundary  layer 
thickness  is  on  the  order  of  i  =  1 .2,  Hence,  the  solutions  are  valid  for  large  values  of  la'.''(2ey))'  -h_,  i.e., 

the  mud  depth  is  large  relative  to  the  wave  boundary'  layer  thickness.  Assuming  that  the  free  surface  is 
uncontaminated,  energy  dissipation  occurs  only  in  the  boundary  layers,  and  is  panitiored  as  follows: 
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where  Uj,  i  =  1.2,3  denote  the  rotational  velocity  components  in  the  direction  of  wave  propagation  in  the  respective 
boundary  layers,  Cd,  is  the  mean  energy  dissipation  rate  in  the  water  boundary  layer  near  the  interface.  Coj  is  the 
mean  energy  dissipation  rate  in  the  mud  boundary  layer  near  the  interface  and  £(>3  is  the  mean  energy  dissipation 
rate  in  the  mud  boundary  layer  near  the  rigid  bottom,  overbar  denotes  time  averaging  and  ”  ±  00 "  implies  "outside" 
the  boundary  layer.  The  various  terms  are  given  by: 

u,=C,exp{-(H-i)tw/(2»'i)}'^(z+h,)}exp(i(kx-at)] 

U2=C2exp{-(l+i)t«/(2i<j)]'^(z+h|)}exp[i(k.x-wti) 

U3=C3exp{-(l  +i)[w/(2v2)]''^(z+h,+h2)}cxp[i(kx  ajt)J, 

C,  =  {-[gka/(wsinhkh,)][(P2/p,)(i'2/v,)''^J/{I+(p,/p,)(vj/»«,)''^l}{{a)-Vgk)-  {(l/2)sinh2kh,J{cothkh, +cothkh,)f(4uVgk)- 
tanhkh,)]}, 

C,=-(P,/P2XV,/V2)‘'’C, 

€3= {-gka/u)Coshkh|  [(wVgk)-tanhkh|]  }/sinhkh2 
and  the  wave  number  is  solved  iteratively  from: 
uVgk  =  [-!-(lMmn)''^J/(2m) 
where 

1  = -(P2/Pi)tanhk(h, +hj)(  I +tanhkh,tanlikh2) 

m = (P2/P 1) + tanhkhitanhkhj 
n = [(P2/P I )- 1  ltanhkh|tanhkh2 


Note  3:  The  rate  of  wave  energy  oissipation  is  computed  by  assuming  mud  as  a  viscous  medium  (Dean  and 
Dalrymple,  1984): 


(A. 5) 


where  the  overbar  denotes  wave-mean  value,  and  h,  and  h,  denote  the  depths  of  wafer  and  mud  layers,  respective!) . 
A  shallow  water  wave-mud  interaction  model  (Jiang  and  Mehta,  1993),  which  describes  a  two-layered  mud/water 
system  forced  by  progressive  surface  waves  of  periodicity  specified  by  wave  frequency,  w.  is  used  to  compute  the 
required  velocity  gradient  inputs. 

Since  the  model  assumes  an  inviscid  upper  water  layer  on  the  premise  that  pressure  and  inertia!  forces  are 
typically  dominant  in  governing  water  motion  and  an  irrotatioiial  flow  field  is  an  apt  description  here,  wave 
dissipation  in  the  water  layer  is  nil.  Furthermore,  by  ignoring  the  velocity  gradient  in  the  vertical  (shallow  water 
assumption),  the  above  equation  simplifies  to: 
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which  may  be  conveniently  divided  into  two  tenns. 


The  appropriate  final  expressions  are  (Feng,  1992): 
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where  A=a/h,,  A^  =  tanh(mhj),  =  cosh(mh2),  =  sinhfmhj),  B,  =  cosh(2mh2),  C,  =  sinh(2mhj), 

hj  =  hj/h, ,  m  =  (-iR^)’'^,  i=(-l)‘^,  and  =  «hj'/«'j  is  the  wave  Reynolds  number,  k  =  kh,  is  the  normalized 
complex  wave  number  given  by: 


i+h2r-[(i*h2r)^-4rh2r]‘'^ 

2rhjr 


(A.IO) 


where  k=k,+ikj,  k,  =  wave  number,  k|  =  wave  attenuation  coefficient,  r=(pj-p,)/p2.  r=  l-tanh(mfiu)/(mh,)  and 
F,=u(h,/g)‘'’  is  the  wave  Froude  number. 


Note  4:  The  solution  is  based  on  a  system  consisting  of  linear  water  waves  propagating  over  a  non-rigid  bed 
characterized  as  a  Voigt  viscoelastic  body  overlain  by  inviscid  water,  with  appropriate  assumptions  to  obtain  an 
explicit  expression  for  wave  attenuation.  By  introducing  a  viscoelastic  parameter,  »>*  =  »>,=»'j-fiG/wp2 .  where  G  is 
the  shear  modulus  of  elasticity,  it  is  shown  that  the  equations  of  motion  for  a  viscoelastic  medium  can  be  reduced 
to  the  form  of  the  linearized  Navier-Stokes  equations  for  a  viscous  fluid.  The  following  assumptions  are  then 
introduced  to  render  the  solutions  explicit: 
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a)  Oscillations  of  the  bed  due  to  the  propagating  waves  are  small  in  comparison  with  the  surface  waves,  as 

in  the  case  of  an  almost  rigid  bed,  or  mathematically,  |  |  ►  I  where  k  is  the  complex  wave  number 

defined  in  Note  3; 

b)  Mud  region  is  assumed  to  be  infinite  in  depth; 

c)  Wave  length  is  assumed  long  compared  to  water  depth. 

The  resulting  expression  for  k,  after  ignoring  higher  order  terms,  becomes; 


<a  ^  P,gG  ip,gi'j 

^gh,  4p2«^{«'2 +GVp2M^)  4pjW(»'2+GVp2W^) 


(A.ll) 


The  mean  rate  of  energy  dissipation  is  then  computed  based  on  the  mean  work  done  on  the  lower  layer  per 
unit  area  of  the  interface  as; 


T 


(A.  12) 


on  z=jj2=0  after  linearizing,  where  T  =  wave  period,  and  p,  and  w,  are  the  pressure  and  vertical  velocity, 
respectively,  in  water.  Upon  time-averaging  and  integration,  the  following  expression  results; 
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is  defined  as  the  amplitude  ratio. 


(A.  13) 


Note  5:  The  dissipation  mechanism  employed  here  is  that  due  to  Yamamoto  and  Takahashi  (l985),  attributed  to 
wave-induced  motion  in  the  soil  bed.  Based  on  a  theoretical  analysis  for  two  elastic  spheres,  Mindlin  and 
Deresfewicz  (1953)  found  that  the  nonlinear  elasticity  and  the  energy  dissipation  of  granular  materials,  e.g.  soils, 
are  due  to  Coulomb  friction  and  are  independent  of  the  forcing  frequency.  In  the  frequency  domain,  the  constitutive 

relation  is  approximately  given  by  the  linear  model,  t  =  G7,  where  r  =  shear  stress,  7  =  shear  strain  and 
G  =  G(1  +i6)  is  the  complex  shear  modulus.  The  real  part,  G,  is  the  dynamic  shear  modulus,  and  5  is  the  specific 
loss  given  by  5=AW/25rW,  in  which  AW  =  energy  loss  per  cycle  and  W  is  the  maximum  elastic  energy  given  by 
W  =  0.507'^,  where  7'  =  shear  strain  amplitude,  both  being  expressed  per  unit  area.  The  effect  of  soil  Coulomb 
friction  on  the  attenuation  of  water  waves  is  approximately  estimated  by  the  analytical  solution  obtained  by 
YamaiiKito  et  al.  (1978)  for  the  quasi-static  response  of  a  homogeneous  half-space  of  porous  elastic  medium  to  water 
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waves.  It  is  assumed  that  the  soil  moduli,  G  and  6,  are  constant  throughout  the  bed.  It  is  noted  that  the  pure  shear 
strain  condition  is  induced  in  the  bed  by  the  water  wave  for  the  entire  phase.  The  magnitudes  of  the  shear  stress, 

r„,  and  the  shear  strain,  y^,  are  only  a  function  of  depth,  z,  and  are  given  by  =  -p^k^ze*^'  and 

7m  =  “(Po/G)koZe'^',  respectively.  Here  po  is  the  amplitude  of  bottom  pressure  given  by  p(,=Piga{l/coshkoh,), 
which  is  the  wave-induced  dynamic  pressure  component  based  on  linear  wave  theory. 

The  energy  loss,  AW,  due  to  soil  Coulomb  friction  is  given  by: 

-o* 

AW  =  2ir5  f  ir  7  dz  (A- 14) 

4  2  " 


which  after  integration  becomes: 


IT  (p,g)^a^5 
2  k^Gcosh^kph, 


(A. 15) 


where  ko  is  the  wave  number  for  a  rigid  bottom.  Accordingly,  the  rate  of  energy  dissipation  becomes: 

^  _  AW  _  (P,g)^  uSa^ 

’  T  4  Gk„cosh’igh, 


Note  6:  Here  the  approximate  analysis  assuming  quasi-static  motion  of  homogeneous  soil  bed  (see  Note  5)  is 
improved  upon  by  utilizing  the  dynamic  theory  for  the  case  of  a  homogeneous  half  space  of  Coulomb-damped 
poroelastic  medium  to  account  for  the  coupling  of  water  waves  and  the  bed  motion  (Yamamoto,  1982;  1983).  While 
two  separate  expressions  for  wave  attenuation  coefficient,  due  to  soil  Coulomb  friction  and  percolation,  respectively , 
are  obtained,  only  the  former  is  used  here  based  on  the  premise  that  typical  permeability  in  mud  is  very  small. 
Substituting  into  the  steady  state  wave  energy  equation  for  a  constant  water  depth,  CjdE/dx  =  -fe,„  and  assuming  an 
exponential  decay  law  then  leads  to: 


1 

4 

GkgCOsh^lCph, 

( I  -  wVwJ  )^  +  5^ 

(A. 17) 


where  wo^=:2Gko^/P2- 

Note  7:  This  semi-empirical  approach  is  based  on  field  measurement  in  East  Bay,  Louisiana  in  water  depth  of  20  m 
where  the  bottom  sediments  are  composed  of  clays  and  silts.  The  instruments  used  included  a  wave  staff,  a  pressure 
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sensor  and  a  bottom-emplaced  accelerometer.  The  general  characteristics  of  the  data  showed  that  the  motions  can 
be  adequately  described  in  terms  of  the  following  harmonic  functions: 

P  =  Po  ^  /o,gAcos(k,x  -  wt)  (A.  18) 

^  ^  MAcos(k,x  -  wt  +  0)  (A.  19) 

where  p  =  wave-induced  bottom  pressure,  =  steady-state  bottom  pressure,  A  =  surface  wave  amplitude,  hjo  = 
depth  of  mud  over  which  motion  occurs,  M  =  proportionality  constant  between  the  amplitudes  of  the  mud  wave 
and  the  pressure  wave;  and  d  =  phase  angle  between  the  crest  of  the  pressure  wave  and  the  crest  of  the  mud  wave. 
The  mean  rate  of  energy  dissipation  is  then  computed  based  on  the  average  energy  transmitted  through  the 
sea/sediment  interface  per  unit  area  over  one  wave  cycle: 

=  I  f  p^dt  (A.20) 

T  dt 

where  dhj  =  an  infinitesimal  increase  in  the  height  of  the  interface.  Using  linear  wave  theory  to  interpret  bottom 
pressures  in  terms  of  surface  wave  height,  and  upon  integration,  the  mean  rate  of  energy  dissipation  due  to  mud 
wave  motion  becomes: 


p.g 

«Msin<f)a^ 

2 

cosh^kjh. 

(A.21) 


where  Values  of  M  and  <t>  can  be  obtained  either  from  a  theoretical  approach  or  field  measurement. 

Note  8:  The  energy  dissipation  per  unit  volume  in  a  linear  viscoelastic  medium  over  one  cycle  is  given  by: 

(tft  «2t  -d’ 

AE  =  f  rydt  =  irylG”  (A. 22) 

..♦s 


where  6'  =  loss  angle,  t  -  shear  stress,  y  =  shear  rate,  y^  =  shear  strain  at  steady  state,  G"  =  loss  modulus. 
Hence,  the  rate  of  energy  dissipation  (per  unit  volume)  becomes: 


2'jr 


AE 


5  2  J  ' 
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where  fi'  =  dynamic  viscosity. 

By  defining  an  average  value  of  y  by  (Schreuder  et  al..  1986): 


{A.23) 
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(A.24) 


<y> 


then, 


jo,' <•,>■>  (A.25) 

The  same  expression  can  also  be  obtained  by  considering  the  response  of  a  viscoelastic  material  to  a  small- 
amplitude  oscillatory  shear.  Given  T=Toe'“,  7=70®  '*“''*’,  G'=T/7=Goe^'=G'-iG"  where  G'  =  complex  shear 
modulus,  G'  =  storage  modulus  and  G"  =  loss  modulus,  the  time-averaged  energy  dissipation  per  unit  volume  over 
a  wave  period,  e,,  reads: 

=  ^  (  Re(r)Re(7)dt  =  (A.26) 

*  ■'o  ^ 

where  Re  denotes  real  part.  By  writing  the  complex  viscosity  as  ,  it  can  be  shown  that  G"-un'  and, 

hence,  the  equivalence  of  the  right-hand  sides  of  Eq.  A.23  and  Eq.  A. 26. 

Assuming  a  uniform  energy  dissipation  per  unit  volume  across  the  depth  of  mud  layer  in  analogy  with  that 
of  uniform  energy  dissipation  in  the  water  column  due  to  breaking-induced  turbulence,  the  mean  rate  of  energy 
dissipation  per  unit  area  can  be  written  as: 
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(A.27) 


Given  h^^h^fx)  as  the  equilibrium  mud  thickness. 
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fin 


(n’  <7>^hj) 


(A.28) 


The  above  expression  can  then  be  evaluated  by  selecting  an  appropriate  viscoelastic  model,  e.g.  Voigt  body,  to 
compute  the  relevant  input  parameters. 
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APPENDIX  B.  A  SUMMARY  OF  SOME  WAVE-MUD  MODELING  STUDIES 


Invesugator(s) 

Water 

Mud 

Modeling 

Brief  Description 

Column 

Depth  Theory 

Gade  (1957, 

1958) 

Inviscid 

Viscous' 

Shallow 

water 

Linear 

Analytical  approach  for  wave 
attenuation;  laboratory  tests  on 
wave  amplitude  decay  and 
interfacial  profiles,  using  a  two- 
fluid  system  (kerosine  overlying 
water-sugar  solution). 

Mallard  and 

Dalrymple 

(1977) 

inviscid 

Elastic' 

Mud 

thickness 

unbounded 

Linear 

Analytical  approach  for  soil 
stresses,  displacement  and  wave 
kinematics.  Soil  inertia  neglected. 
No  energy  di.s.sipation  since 
viscous  property  of  mud  not 
considered. 

Dalrymple  and 

Liu  (1978) 

Viscous 

Viscous' 

Linear 

Analytical  approach  for  wave 
attenuation.  Model  valid  for  any 
water  depth  and  both  deep  and 
shallow  mud  layers;  boundary- 
layer  approximation;  explicit 
solutions. 

Dawson  (1978) 

inviscid 

Elastic' 

Mud 

thickness 

unbounded 

Linear 

Same  as  Mallard  and  Dalrymple 
(1978),  but  soil  inertia  taken  into 
account. 

Hsiao  and 
Shemdin  (1980) 

Inviscid 

Viscoelastic' 

(Voigt) 

Linear 

Analytical  approach  for  wave 
attenuation. 

MacPherson 

(1980) 

Inviscid 

Viscoelastic' 

(Voigt) 

Linear 

Analytical  approach  for  wave 
attenuation. 

Yamamoto 
(1982,  1983); 
Yamamoto  and 
Takahashi 
(1985) 

Inviscid 

poroelastic” 

Multi¬ 

layered 

mud 

Linear 

Based  on  nonlinear  poroelastic 
theory  that  includes  nonlinear  and 
imperfect  elasticily  of  soil 
skeleton.  Bed  respon.se  nonlinear 
and  dynamically  amplified; 
primary  dissipation  mechanism 

internal  Coulomb  friction  between 
grains.  Arbitrary  distribution  of 
continuou.siy  varying  soil 
properties  with  depth.  Qua.si-static 
response  more  appropriate  for 
hard  bed  .such  as  sands,  while 
dynamic  response  for  soft  beds 
such  as  clays. 


Mehta  and  Maa 

Viscous 

Viscoelastic' 

Multi¬ 

Linear 

Analytical  approach  with 

(1986);  Maa 

(Voigt) 

layered 

laboratory  tests  on  wave 

and  Mehta 
(1987,  1989) 

mud 

attenuation  and  mud  motion. 

76 


Investigator(s) 

Water 

Column 

Mud 

Modeling 

Depth 

Theory 

Brief  De.scription 

Shibayama  et 
al.  (1986) 

Viscous 

Viscous' 

Linear 

Analytical  approach  using  explicit 
velocity  field  of  Dalrymple  and 

Liu  (1978)  and  laboratory  tests  on 
mud  mass  transport  and 
attenuation  rate. 

Suhayda  (1986) 

Inviscid 

Viscoelastic' 

(Nonlinear) 

Mud 

thickness 

unbounded 

Linear 

Analytical  approach  based  on  the 
explicit  results  of  MacPherson 
(1980);  assumptions  of  small 
interfacial  oscillations  and  infinite 
mud  thickness.  Nonlinear 
expressions  relating  shear  strain  to 
shear  stress  and  damping  ratio. 

Mei  and  Liu 
(1987) 

Inviscid 

Viscoplastic' 
(Bingham, 
thin  layer) 

Shallow 
water; 
thin  mud 
layer 

Linear 

and 

Solitary 

Analytical  approach  for  wave 
attenuation.  Laminar  constitutive 
model  and  interfacial  shear  stress 
<  <  mud  yield  stress.  Predicts 
continuous  to  intermittent  mud 
motion. 

Tsuruya  et  at. 
(1987) 

Viscous 

Viscoplastic' 

(Bingham) 

Multi¬ 

layered 

mud 

Linear 

Analytical  approach  by  extending 
Dalrymple  and  Liu’s  (1978) 
viscous-fluid  model  to  Bingham 

Foda 

(1989) 


Inviscid 


Stratified 

Elastic' 


Deep  water; 
mud 

thickness 

unbounded 


Liu  and  Mei 
(1989) 


Viscous 


Viscoplastic' 

(Bingham) 


Shallow 
water; 
thin  mud 
layer 


Solitary 


fluid.  Laboratory  tests  on  mud 
mass  transport  and  wave 
attenuation.  Equivalent  viscosity 
that  takes  into  account  the  non- 
Newtonian  mud  rheology  is 
introduced  in  an  approximate 
fashion  into  the  Newtonian 
equations  of  motion  and  solved 
iteratively. 

Analytical  approach  for  wave 
attenuation.  Wave  energy  lost  to 
highly  dissipative^  sideband 
oscillations,  which  in  turn  lose 
energy  due  to  viscous  dissipation 
whose  action  is  enhanced  by  the 
interaction  of  standing  elastic 
shear  waves  with  viscous 
boundary  layer. 

Analytical  approach  that  includes 
interfacial  shear  through  boundary 
layer  approximation.  Wave 
dissipation  due  to  bottom  mud 
shear  layer  and  turbulent 
interfacial  stress.  Predicts 
continuous  to  intermittent  mud 
motion. 


77 


Investigator(s) 

Water 

Column 

Mud 

Modeling 

Depth 

Theory 

Brief  Description 

Jiang  and 

Zhao  (1989) 

Viscous 

Viscous* 

Shallow 

water 

Solitary 

Analytical  approach  and  laboratory 
tests  on  wave  attenuation  based  on 
boundary  layer  approximation. 

Chou  (1989) 

Viscous 

Viscous/ 

Viscoelastic/ 

Elastic* 

Three¬ 

layered 

mud 

Linear 

Analytical  approach  for  wave 
attenuation  and  fluidization.  A 
generalized  viscoelastic  model 
where  the  bed  consists  of  two 
layers  with  viscous/viscoelastic 
properties  overlying  a  half-plane 
of  viscoelastic/elastic  muds,  and  a 
viscoelastic  boundary  layer 
approximation. 

Shibayama  et 
al.  (1989) 

Viscous 

Viscoelastic* 

Multi¬ 

layered 

mud 

Linear 

Analytical  approach  and  laboratory 
tests  on  mud  mass  transport. 

Sakakiyama  and 
Bijker  (1989) 

Viscous 

Viscous* 

Linear 

Analytical  approach  based  on 
Dalrymple  and  Liu’s  (1978) 
viscous  fluid  model  with 
laboratory  tests  on  mud  mass 
transport. 

Jiang  et  al. 

(1990) 

Viscous 

Viscous* 

Shallow  water 

Cnoidal 

(1« 

order) 

Analytical  approach  for  wave 
attenuation  based  on  boundary 
layer  approximation.  A  more 
general  model  that  includes  those 
of  Jiang  and  Zhao  (1989)  for 
solitary  waves  and  Dalrymple  and 
Liu’s  (1978)  for  their  explicit 
solutions  based  on  boundary  layer 
approximation  under  shallow 
water  waves. 

Shibayama  et 
al.  (1990) 

Viscous 

Viscoelastic 

and/or 

viscoplastic* 

Linear 

Extension  of  the  model  of 

Shibayama  et  al.  ^989)  to  include 
viscoplastic  bed  when  the  shear 
stress  is  greater  than  the  Bingham 
yield  .stress. 

Feng  (1992) 

inviscid 

Power 

Fluid* 

Shallow  water 

Linear 

Analytical  approach  and  laboratory 
tests  on  wave  attenuation  and  mud 
fluidization.  Power  fluid  (Sisko) 
model  obtained  empirically. 

Isobe  et  al. 

(1992) 

Inviscid 

Non- 

Newtonian* 

Linear 

Numerical  approach  for  wave 
attenuation  and  mud  mass 
transport  in  which  an  empirical 
rheologic  mud  model  based  on 
multi-linear  approximation  to 
experimental  stress-strain  rate 
curves  is  derived. 
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Investigator(s) 

Water 

Column 

Mud 

Modeling 

Depth  Theory 

Brief  Description 

Jiang  (1993) 

Viscous 

Viscoelastic' 

Non¬ 

linear 

Analytical  approach  and  laboratory 
tests  on  wave  attenuation  and  mud 
mass  transport. 

Cueva  (1993) 

Viscous 

Viscoelastic 

(Voigt) 

Linear 

Analytical  approach  for  wave 
attenuation  based  on  boundary 
layer  approximation 

*  Mud  is  assumed  as  a  single-phase  continuum. 

Two-phased  system  consisting  of  soil  skeleton  and  pore  fluid  in  relative  motion. 
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presented.  Models  for  water  waves  over  mud  are  first  presented,  then  compared  to  laboratory  data.  Field  data  are  used  to 
make  example  applications  of  the  models. 
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